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In group-testing, a set of x units is taken from a total starting set of N 
units, and the x units (1 ^ x ^ N) are tested simultaneously as a group 
with one of two possible outcomes: either all x units are good or at least one 
defective unit is present (we don't know how many or which ones). Under 
this type of testing, the problem is to find the best integer x for the first test 
and, to find a rule for choosing the best subsequent test-groups (which may 
depend on results already observed), in order to minimize the expected total 
number of group-tests required to classify each of the N units as good or 
defective. It is assumed, that the N units can be treated like independent 
binomial chance variables with a common, known probability p of any one 
being defective; the case of unknown p and several generalizations of the 
problem are also considered. 

I. SUMMAKY 

A finite number, N, of units are to be tested in groi ps. A "group-test" 
is a simultaneous test on x units (x to be chosen so that 1 ^ x ^ N) 
with only two possible outcomes: "success," indicating that all x units 
are good; and "failure," indicating that at least one of the x units is 
defective (we don't know how many or which ones). The problem is to 
define a simple and efficient procedure (or an optimal procedure) for 
separating all the defective units from the good units — efficiency being 
defined in the sense of minimizing the expected number of group-tests 
required. Each unit is assumed to represent an independent observation 
from a binomial population with a common known a priori probability, 
q, of being good and p = 1 — q of being defective. (The case of q un- 
known is briefly treated in Section X.) 

A procedure (or decision rule), Ri , which describes a mode of action 
for any given value of q, is proposed and compared with several other 
procedures applicable to the same problem. The procedure R\ is simple 
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in the sense that at any time, t, the experimenter must separate the 
units not yet proven to be good or defective into only two sets; units 
within either of these two sets need not be distinguishable. If it is given 
that the identification of units within the group being tested is econom- 
ically impractical or impossible, then the procedure Ri is conjectured to 
be optimal for all values of q. 

Explicit instructions for carrying out Ri are given for iV = 1(1)16 for 
all q and for N = 17(1)100 for the particular values q = 0.90, 0.95 and 
0.99. Exact formulae for the expected number of group-tests required 
under R { are given in Table IV B for all q and for N = 1(1)12; numerical 
results for q = 0.90, 0.95 and 0.99 are given in Tables V A, V B and V C 
for AT = 1 ( 1 ) 100. Other numerical comparisons are made in Tables II A 
and II B. [Tables II through VIII appear at the end of this paper.] 

Another procedure, R 2 , which is simpler to compute and compares 
favorably with Ri , is defined in Appendix A in terms of information 
theory concepts. 

Several different directions for generalization of the problem and cor- 
responding generalizations of the procedure Ri are considered in Section 
XI. Industrial applications are mentioned, in addition to the known 
application to blood testing. 

II. INTRODUCTION 

A problem which has hitherto been considered only in connection 
with blood-testing applications 1 ,2,3 can be shown to have industrial 
applications, and these have focused interest on a more general treat- 
ment of the problem. During World War II, a great saving was accom- 
plished in the field of blood testing by pooling a fixed number of blood 
samples and testing the pooled sample for some particular disease. If the 
disease was not present, then several people were passed by a single test; 
if the disease was present, then there was enough blood remaining in 
each blood sample to test each one separately. The amount of time, 
money and effort saved by such a procedure depends on how rare the 
disease is in the population of people being tested. In this application, 
the total number of people to be tested was regarded as unknown and 
very large. 

The goal of the problem treated here is the same — namely, to separate 
the defective units from the good units with a minimal (or approxi- 
mately minimal) number of group-tests. This problem differs from the 
blood-testing problem in the following respects: 

i. The population size N (number of people to be tested) is known 
at the outset. 
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ii. The number of units in each group-test (pooled blood sample) 
is not necessarily constant. 

iii. If a group-test fails (the disease is present) we do not necessarily 
test each item separately. 

In practice, the simplicity of the procedure deserves some considera- 
tion. The proposed procedure /?i defined in Section III, after having 
been computed and described explicitly in advance of any experimenta- 
tion, is in some sense no more complicated than the blood-testing pro- 
cedure described above ; this is explained in Section V. 

Some typical industrial applications are: 

1. It is desired to remove all "leakers" from a set of N devices. One 
chemical apparatus is available and the devices are tested by putting 
x of them (where 1 ^ x ^ N) in a bell jar and testing whether any of 
the gas used in constructing the devices has leaked out into the bell jar. 
It is assumed that the presence of gas in the bell jar indicates only that 
there is at least one leaker and that the amount of gas gives no indica- 
tion of the number of leakers. The a 'priori probability, q, of a unit being 
good is given by the records of similar units tested in the past. 

2. Paper capacitors are tested at most n at a time, and each test in- 
dicates by the presence or absence of a current whether or not there is 
at least one defective present. For given n and q and given cost of unit 
manufacture, should the operator throw away a whole set of n units if 
it contains at least one defective? If not, how should he proceed to sort 
out the defective units to minimize the expected number of tests re- 
quired? If the cost of a group-test and the cost of producing a unit are 
known, a related problem is to find a procedure which minimizes the 
total cost (including testing costs) of producing a good unit. 

3. Christmas tree lighting problem. A batch of n light bulbs is electric- 
ally arranged in series and tested by appljdng a voltage across the whole 
batch or any subset thereof. If this is to be done on a routine basis, what 
procedure should be used to minimize the expected number of tests re- 
quired to remove all the defective light bulbs, assuming the value of q 
is given? 

4. A test indicates whether or not there is at least one good unit present 
in a batch of n, without indicating which ones or how many are good. 
Given q, what procedure should be used to remove the good units? This 
dual problem, which is useful in salvaging good components on a routine 
basis, is mathematically equivalent to those above, if the definitions of 
good and defective are interchanged. 

A procedure R t is defined to solve the above problems, and is compared 
with several other procedures for the same problem. A procedure R> , 
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based on maximizing the information in each group-test, is defined in 
Appendix A. Another procedure, R 3 , which does not allow any recom- 
bination, is defined in Appendix B. Two "halving procedures," which 
can be carried out without knowing the true value of q, are defined in 
Appendix C. Procedures R 7 and R 6 are the best procedures that can be 
obtained by the methods of Dorfman 1 and Sterrett, 3 respectively, for 
small population sizes. For N = 4, 8 and 12 and various q values, Table 
II A gives a numerical comparison of the expected number of group-tests 
required for all these procedures. 

Different directions of generalization, some of which are discussed in 
Section XI, are the following: 

1. Two (or more) different kinds of units with (say) known proba- 
bilities q y , 52 of a unit being good are present, and the two different kinds 
can be put into the same test group. 

2. Two (or more) experimenters may be working on a single set of N 
units by carrying out simultaneous, parallel group-tests and cooperating 
in such a way as to minimize the time required to accomplish the task. 

3. The restriction is sometimes applied (particularly in blood-testing) 
that any one unit can be included in at most k group-tests; here the goal 
is to minimize the expected number of group-tests subject to this re- 
striction. For fc = 2 the proposed procedure is necessarily based on the 
method of Dorfman; 1 i.e., if a group-test fails, then the units therein 
are all tested individually. 

4. Various generalizations appear if it is assumed that each test on x 
units gives three (or more) different possible results. For example, a test 
could indicate that either (a) all are good or (b) all are defective or (c) 
there are at least one good unit and at least one defective present. 

5. A unit can be defective in either of two ways (e.g., electrical or 
mechanical) with the two a priori probabilities of being defective as- 
sumed to be independent but not necessarily equal. If there are two 
different tests corresponding to the two types of defectives, then, in addi- 
tion to deciding the next test-group size, it may be necessary to decide 
which test to use next. 

G. For positive continuous chance variables with a known distri- 
bution (like weight) the following problem is analagous. It is desired to 
separate N units into two groups according as the weight per unit is 
less than or greater than a constant (say, unity). Any number of units 
can be included in a single weighing. The problem is to accomplish the 
separation in a minimal number of weighings, assuming that the indi- 
vidual weights are independent observations from the common known 
distribution. 
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Many of these generalizations will be omitted from this paper and 
treated separately. 

III. THE PROCEDURE Ri 

The procedure Ri is denned implicitly by a pair of recursion formulae 
and boundary conditions, but first we shall need some definitions and 
preliminary results. The units proven to be good and the units proven 
to be defective are never used in subsequent tests. Aside from such units, 
this procedure requires that at every stage the remaining units be sep- 
arated into at most two sets. For one set of size m ^ 0, which we call 
the defective set, it is known that it contains at least one defective unit; 
for the other set of size n — m ^ 0, which we call the binomial set, our 
a posteriori knowledge is, so to speak, in the original binomial state; i.e., 
given the past history of testing, the units in the binomial set act like 
independent binomial chance variables with a common probability p of 
being defective. For the defective set, the conditional probability that 
Y, the number of defectives present, equals y is 



Ch 



q m 

Pr [Y = y\Y £ 1} ^ x y_ - (y = 1,2, • • -,wi). K) 

If X denotes the number of defectives present in a subset of size x ran- 
domly chosen from the defective set, then 

™\„v™-v( m - v 
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Pr{Z = 0|F^lJ=E^ VW^ gd-g ) . (2) 

1 ' £i 1 - q m /m\ 1 - q m 

Before defining the procedure it is convenient to prove a lemma in a 
more general setting. Let T(r,) (i = 1, 2, • • • , t) denote a test on n 
units (1 ^ ri ^ n) such that there are only two mutually exclusive 
possible outcomes: a "failure," indicating that there are at reast r,- de- 
fectives present, and a "success," indicating that at most r, — 1 of the 
units in the test are defective. In Lemma 1 we consider any integers r,- , 
ro with 1 ^ r, ^ n(i = 0, ], 2, • ■ ■ , t), but the most important appli- 
cation is the case n — r t = ■ • • = r t = 1. Let & be any set of units, 
and let (B,(i = 1,2, • • • , t) denote sets not necessarily disjoint from 
one another but such that each is disjoint from Ct; the case t = 1 is 
the one used for procedure R\ . At the outset, all units are independ- 
ently and binomially distributed with a common probability p of be- 
ing defective. 
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Lemma 1: If a test T(r { ) on a + (B* produces a failure for (i = 1, 
2, ■ • • , t) and another test T(r ) on a also produces a failure, then for 
r ^ max r t - the conditional distribution associated with all the units in 
the sets <B<(» = 1, 2, • • • , t), given both conditions above, is exactly 
the same as the original binomial distribution. 

Proof: Let A and Bi denote the chance number of defectives present 
in a and CB,- , respectively. For the jth set (By the conditional probability 
P of interest is 

P = Pt{Bj ^ b\A + B;^ n(i = l,2,---,t),A ^ r ). (3) 

Since r Q ^ max n(i = 1, 2, ■ ■ ■ , I), the condition A ^ r implies that 
A -+- #»• ^ r,- , and hence 

P = Pr{B/ ^ 6|i4 ^ ro}. (4) 

Since <B; and a are disjoint, it follows that Bj and A are independent 
and hence, from (4), 

P = Pr{Bj^b}, (5) 

which proves the lemma. 

Let Gi(m, n; q) = Gi(m, n) denote the expected number of group- 
tests remaining to be performed if the defective set is presently of size 
m, the binomial set is presently of size n — m, the a priori probability of 
a good unit is the known constant q and the procedure Ri is used. For 
the special case m = we use the symbol Hi(n; q) = Hi(n). The values 
of m and n vary as the procedure is carried out; at the outset, m = 
and n = N. It will also be convenient to refer to the G-situation or 
G(m, n) -situation if m > 1 and to the H -situation or H(n) -situation if 
m = 0. 

Recursion Formulae Defining Procedure R\ 

If x denotes the size of the very next group-test, then we write for 
any situation with m = 

H x { n ) = 1 + min {q'H^n - x) + (1 - q x ) Gi(x, n)\, (6) 



l<x<n 



and, with the help of (2) and Lemma 1, we write f or n ^ m ^ 2 
Gi(m,n) = 1 + min 1( ? ~ \ ) <h(m - x,n - x) 



(7) 
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The boundary conditions state that for all q 

//,(()) = and Gi(\,n) = H 1 (n - 1 ) for n = 1,2, ■• • . (8) 

In (G) and (7) the constant 1 represents the very next group-test of 
size x and the expression in braces is the conditional expected number 
of additional group-tests given x. It follows from (6) and (8) that 
//,(1) = 1 for all q. 

Remark 1 : To justify writing Gi(z, n) in (7) we make use of Lemma 1 
with t = 2, r 2 = 0, i\ = ra = 1 ; we take CI + (Bi as the defective set of 
size m and ft as a subset of size x < m. Then, by Lemma 1, if the sub- 
set ft of size x is shown to contain at least one defective, the a posteriori 
distribution associated with the m — x units in Bi is exactly binomial. 
These are then mixed or recombined with the n — m "binomial units," 
giving a total of n — x binomial units, and this justifies the expression 
Oiix, n) in (7). 

Remark 2: These two recursion formulae, together with the boundary 
conditions, allow one to compute successively for any q the functions 
Oi(2, 2), Hi(2), (7x(2, 3), (?i(3, 3), 7^(3), G l (2, 4), G 1 (3, 4), d(4, 4), 
//i(4), • • • to any desired value of m and n. 

Remark 3: The integer x which accomplishes the minimization in (6) 
and (7) for each situation characterized by the integers m and n is par- 
ticularly important, since this is the size of the next test to be run ac- 
cording to the procedure R\ . These integers x = x H (n; q) and x = x {m, 
n; q) implicitly define the procedure Ri . An illustration of how the pro- 
cedure Ri is to be carried out is given in Section IV. 

Remark 4: If in > 1, then it is assumed in (7) that a subset of size x 
with 1 ^ x < m will be taken from the defective set without mixing it 
with units from the binomial set. It follows from (G), (7) and (8) that 
any lack of optimality can only arise from this "no mixing" assumption. 
This assumption was used in the derivation of the algorithm (7) (see 
Remark 1 above). It will be noted in Section XIII that, when all the 
units are individually identified, then by dropping this assumption, an 
improvement to the procedure R\ for high values of q can be found. A 
specific example of a modification and improvement of the procedure 7?i , 
which drops the "no mixing" assumption at the expense of more com- 
plication will be more thoroughly discussed in a separate paper. 

IV. ILLUSTRATION OF THE PROCEDURE R\ 

Suppose we start with .V = 12 units and it is given that q = 0.98. 
Referring to the column headed //i (12) in Fig. 3, we find that the first 
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test-group is of size x = 12; i.e., we start by testing all 12 units. If a 
success occurs, the experiment is over; if a failure occurs, then accord- 
ing to the column headed Gi(12, n) of Fig. 4 the next test group is of 
size x = 4 chosen at random from the 12. Similarly, we continue along 
one of the sample paths shown in Fig. 1. The complete "tree" is not 
shown here, but testing continues in a similar manner and the specific 
details can be obtained from Figs. 3 and 4, which appear at the end of 
this paper. 




Fig. 1— Initial part of the tree for procedure R t for q = 0.98 and starting with 
an H{\ 2) -situation. 

It is obvious that the above procedure terminates in a finite number of 
steps. In fact, it can be shown for procedure Ri (proof is omitted) that 
the maximum number M(n) [M(m,n)] for any .//-situation [(/-situation] 
occurs when q is close to unity and the n unanalyzed units are all defec- 
tive. It follows easily that 

M(n) = (n + 1) [1 + «(»)] + 1 - 2 l+aM , 

M(m,n) = a{m) + n [1 + a{n - 1)] + 1 - 2 1+a(n_1) , (m > 1) 
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where a(z), for any positive integer z, is denned by 



!«(:) 



^ z < 2 H * W . 



(ID 



For the example above, a( 12) = 3 and Af (12) = 37 ; it is interesting to note 
that the constant term in Hi(n) expressed in powers of q in Table IV B 
is also the length of the longest "chain" (that is, 37 for n = 12) in the 
"tree" used for the interval of ^-values ending at unity. Although the 
maximum number is so large, the expected number of tests is only 2.07. 
This is explained partly by the fact that the procedure terminates after 
one group-test with probability 0.7847. A table of such probabilities for 
the number of tests T required when q = 0.98 and N = 12 is given 
below : 



T 


l 


2 


3 


4 


5 


6 


7 


8 


9 


10-37 


Proba- 
bility 


0.7847 











0.0801 


0.1124 


0.0003 


0.0016 


0.0134 


0.0075 



If we assign the probability 0.0075 to T = 10, we obtain an estimate of 
//i(12) (namely, 2.070) which is a lower bound. The exact value, 2.073, 
can be obtained from the formula for Hi(l2) in Table IV B. Similarly, 
an estimate of the standard deviation is computed to be a = 2.1, and 
this also is easily shown to be a lower bound. 

It is interesting to note that, if the starting number N is exactly a 
power of 2 and q is large, the procedure Ri starts off the same as a 
"halving" procedure. Such a procedure R\ is defined in Appendix C for 
any N, and it has the property that it can be carried out without know- 
ing the true value of q. To compare the means and standard deviations 
of Ri and Ri , we consider the case AT = 6, where M = 14 for both Ri 
and Ri . For any q > 0.844, the expectation under Rx can be put in the 
form 

EiT-M = l(g 6 ) + 4(3pr/) + 5(:W + 2pV) + 6(2pV) 
+ 7(7pV + 2pV/) + 8(4pV + 5pV) + 9(8pV + 3pV) 
+ 10(5pV + 5pV) + H(3pV + 2p & q) 
+ 12(3pV + V\) + 13(pV + 2p 5 ?) + U(p b q + p 6 ), 

where, for each term, the expression in parentheses is the probability 
that T takes on the value of the associated integer coefficient. For any q 
the corresponding expression under Ri is 



(13) 
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E(T;R<) = l(r/) 4- 4(3pg 6 ) + 5(3p? 5 + 2p 2 g 4 ) + 6(pV + pV) 
+ 7(7pV) + 8(5pV + 5pV) + 9(6pV + 2pV) 

+ io(7pV + 4pV) + h(pV + 4p¥ + p 5 q) 

+ 12(4pV + 2p 5 ?) + 13(pV + 2p 5 ?) + 14(p 5 ? 4- p 6 ). 

The procedure Ri is better than R A (at least for </ > 0.844) since 

E(T;R A ) - EiT-M = 2p 2 q + 5p s q 4- 4pV + p 5 q ^ 0. 

[The fact that each term in (14) is positive indicates that Ri is better 
than Ri , even if we know how many defectives are actually present 
among the N units.] Some numerical comparisons of means and standard 
deviations are given in Table I. The maximum difference for N = 6 be- 
tween E(T; Ri) and E(T; R t ) occurs at q = 0.844, and is equal to 0.0379. 
For N = 6, the procedure Ri appears to have a variance smaller than 



(14) 



Table I — 


Comparison of Procedures Ri 


\nd Ri for N = 6 




q = 0.85 


g = 0.90 


q = 0.95 


q = 0.98 


q = 0.99 


E(T; R{) 
E(T; Ii t ) 


3.8118 

3.8472 


2.9434 
2.9605 


2.0092 
2.0136 


1.4133 

1.4141 


1.2083 

1.2085 


a(T; R^ 
a(T; R t ) 


2.536 
2.593 


2.269 
2.304 


1.762 
1.776 


1.179 
1.183 


0.850 
0.851 



that of 7^4 for all q < 1. A more complete comparison of E(T\ Rx) with 
that of several other procedures is given in Table II A. 

V. THE SIMPLICITY OF R\ 

It will be shown in this section that, for any given q and any situation 
G(m,n), the appropriate x [i.e., the integer which accomplishes the mini- 
mization in (7)] does not depend on n. A somewhat simpler method of 
computing x is given and a new function of m alone is introduced to re- 
place Gi(m,n) in the definition of the procedure R] . For any m and any 
pair of integers (.r, x + 1) both possible under Ri , there is always a 
single dividing point qa(x) = q G (x, x -{- 1; m) that separates the interval 
for x from the interval for x + 1. (This property was observed for 
m ^ n ^ 16 and is treated as a conjecture for all m and n in Section VII.) 

According to Remark 4 in Section III, the procedure Ri for m > 1 is 
to "break down" the defective set. This "breaking down" is continued 
until a single unit is established to be defective and removed. Instead 
of randomizing the order of the units in the defective group again and 
again before each test group is selected, it will be convenient to assume, 
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without affecting the properties of the procedure /?i , that the order is 
randomized only once at the outset. t Units or groups of units removed 
later are then to be taken in that order. If the ith unit (in that order) is 
the first defective unit, then the "breaking down" mentioned above 
leads to an //-situation with n — i binomial units, and the converse also 
holds true.} 

It is convenient to introduce Fi(m, q) = Fi(m) defined as the expected 
number of group-tests required to "break down" a defective set of size 
m and for the first time reach an //-situation when q is given and the 
procedure R t is used. Then Fi(m) clearly does not depend on n and the 
above argument permits us to write 

Gi(m,n) = F 1 (m) + f=-^— ) Z <T'//i(n - 0- (15) 

\1 — q m / ,=i 

For algebraic simplicity we let 

Ch*{m, n) = (p^-) Gx (m, n) and Fi*(m) = (^r") Mm). (10) 
Then (7) and (15) take on the simpler forms 
Gi*(m, n) = 2 q'~* 



17) 



+ min {qGi*{m - x,n - x) + Gi*(x,n)\, 

1 <i<m— 1 



G,*(w,n) = F x *(m) + £ q^fhin - i). (18) 

1=1 

Substituting ( 18) in (17), the three summations cancel and the result is 

I'\*(m) - Ef + min Jr/F,*(w - x) + Fi*(x)} 3 (19) 

which does not depend on n. The boundary condition, Fi*( 1 ) = for all q, 
also does not depend on n. It is clear from this derivation that (19), 
which does not depend on n, must define the same integer values 
x = x (?n;q) as (17) or (7). This proves the following theorem. 

t It should be pointed out that even this single randomization at the outset 
can be disregarded in carrying out the procedure R\ if there is no doubt about the 
assumption of independent chance variables or if the units are already well-mixed 
in the process of delivery to the experimenter. 

X It follows from the above that, for any procedure which "breaks down" the 
defective set in the above-mentioned manner (including a method of testing units 
from the defective set one at a time until a defective unit is found), the expected 
number of good units eliminated between a G(nt,n) -situation and the next //-situa- 
tion is q/p — >nq'"/ (I — </'"), and the number of defective units eliminated is 
always exactly one. 
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Theorem 1 : For any (7-situation with n ^ m > 1 and any q, the size 
of the next test group, denned implicitly by (7), does not depend on n. 

This result simplifies the explicit instructions needed to describe the 
procedure. Thus the two diagrams, Figs. 3 and 4, describe the procedure 
Ri for all values of q and for any N ^ 16. 

Equations (15) and (16) can also be substituted in (6), yielding 



fain) = 1 + min \q x Hi(n - x) 

+ (1 - q) [fx*(x) + g t'H.in - *)]) , 



which, together with (19), gives a pair of "one-dimensional" recursion 
formulae for denning Ri instead of the "two-dimensional" set, (6) and 

(7). 

Remark 5: It .should be pointed that, if one were to ask for a procedure 
that "breaks down" the defective set in as small an expected number of 
group tests as possible, then one would write (19) as one of the basic 
recursion formulae defining the procedure. This shows that Ri "breaks 
down" the defective set and returns to an //-situation in a minimal num- 
ber of tests. 

VI. SOME PROPERTIES OF Ri FOR q CLOSE TO UNITY 

For any G-situation with m > 1, consider the effect of increasing q. 
It is easy to see by an induction argument that G\(m, n) is a strictly 
decreasing function of q. The function H^n) is also strictly decreasing 
unless the value of q is such that the procedure Ri tests all units one at 
a time, in which case Hi(n) is constant. The "tree" remains the same 
in an interval as q increases, and changes only when it becomes more 
efficient under Ri to increase the size of some test group in the "tree"; 
i.e., if we proceed down the "tree" along any path, the first change en- 
countered, if any, will be an increase in some test group size. It there- 
fore seems reasonable to expect (in both G- and //-situations) that the 
largest value of x assigned by R t (say, x max ) occurs in an interval of q- 
values ending at unity. This unproved assertion that under Ri large 
values of x are associated with large values of q is an immediate conse- 
quence of Conjectures 1 and 2 stated in Section 7. 

For fixed m > 1, let the integers a(m) and (i(m) 

[a(m) ^1, ^ j8(OT) < 2 o(m) ] 
be defined by 

m = 2 a(m) +/3(m). (20a) 



(21) 
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Under the assumption that R t assigns x max . in an interval of q values end- 
ing at unity, it will be shown for any G-situation that in this interval 
the value of x = x max is given by 

| 2 «(»)-i for 2 «eo ^ m g 3.2-c-)- 1 

' r "' ax (m - 2 a(m) for 3-2 a(M - 1] ^ m < 2 B( "° +1 . 
As a corollary it then follows that, under #1 (G-situation), 

1 S x„« § » (22) 

Also in the above-mentioned interval we have 

ljZXj + ff "-*<*) 1^-f—J. (23) 

Proof of {21), (22) and (23): In (19), since z ^ m - 1 and 
m — # ^ m — 1, we can use the induction hypothesis (23) to obtain 

m f /l — a"' -I \ 
Fi*(«l) - 2 fl' -1 = min sg x a(wi - x) ( — ) 

i-1 ig*g»-i [ \ 1 — <? / 

+ ,-„,„., (L^r;') + .,,, (L^i) (21, 

For q close to unity, the right side of (24) is equivalent to minimizing 

Q(x) = [xa(x) + 20(a)] 4- [(m - rr>(m - x) + 2/3 (m - *)]. (25) 

For the moment, let us utilize the symmetry of Q(x) for 2 ^ x ^ ?n — 2 
and limit our considerations to x ^ m/2 and ra ^ 4. 

Let [rn/2] denote the largest integer less than or equal to m/2. For 
x = [m/2], [m/2] - 1, • • • , we consider Q(x + 1) - Q(x) and distin- 
guish several cases according as 

a. a(x — 1) = a(x) = a(m — x) ^ a(m — x + 1); 

b. a(;C - 1) 4- 1 = a(x) = a(m - x) ^ a(m — a; + 1); 

c. a(x — 1) = a(z) = a(w — X — 1) 

= a(w - x) - 1 = a(m - X + 1) - 1; (2G) 

d. a.(.r —1)4-1 = «(•?) = a(wi — x — I) 

= a(m — x) — 1 = a(w — x + 1) — 1; 

e. m ^ 2x- 4- 1 and a(x) < a(m — x — 1). 
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Using the fact that, for any integer y > 1, 

y[a(y) - a(y - 1)] + 2\fi(y) - fry - 1)] = 2, (27) 

we obtain for all the cases in (26) the result, 

Q(x + 1) - Q(x) = a(x) - a(m - x) ^ 0. (28) 

In particular, the value is zero only when (26a) holds and a(x) = 
a(m — x). Hence, Q(x) is a nonincreasing function of the integer x 
for 1 ^ x ^ w/2 and is constant for x ^ x ^ m — Xq , where Xq is de- 
fined by 

a(x Q — 1) + 1 = a(xfi) = a(m — x — 1) 

or (29) 

a(x ) = a(m — .To — 1) = oc(m — Xo) — 1, 

whichever gives the maximum. These imply that 

Xo = 2 «(«o) ^ m _ Xo _ j < 2 «t*o)+i = 2 . r „ 

or (30) 

.To < TO - .To = 2 a(XD)+I ^ 2To . 

Both lead to the same results; namely, 



It follows that 



.5- ^ -To < ^ . (31) 



to , . 2to 

3 s x '"« s T 



and, from (29), for any integer x with .t % x ^ to — .t , the values of 

a(.r) and a(m — .t) can differ by at most unity. 

For any x with .t ^ a; ^ to — .t we consider two cases according as 
I a(x) — a(m — x) I is zero or unity. Let the expression in braces in (24) 
be denoted Ci*(x). For integers x with#o ^ x S m — x the number of 
terms in Ci*(x) is constant and the expression C*(x) is to be minimized 
by making more powers of q larger. 

Case 1. We can write 

[1 »'"1 [~1 20(m-x)-l 



(32) 

T .„_, II - - n ipw I 



L I"? J' 
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and the problem is to find the x which minimizes Ci*(x) for large values 
of q (i.e., to find .r„ iax ). Since a(x) = a(m — x), (3(m — x) — /3(.r) = 
m — 2.r and, using the fact that x max ^ ra/3, we consider only integers 
x ^ ra/3 and obtain 

m - x ^ 2(ra - 2x) = 2[&{m - x) - P(x)], (33) 

and hence 

m - 2j8(m - .r) ^ .r - 2/3(.r). (34) 

It follows that the last term in (32) has lower powers of q than the 
previous term and that (32) is minimized by setting 0(x) = 0, i.e., by 
taking x max to be a power of 2. To complete the proof of (21) and (23), 
we note that 

m = x + (m - x) = 2 a<r) + 2 a( '"~ l) + fi(m - x) 

= 2 a <* )+1 + fi(m - x) t { ) 

so that a(m) = a(x) + 1, /3(ra) = /3(ra — x) and, since .r = .r Minx is a 
power of 2, we have 

x max = 2 a(m) -\ (36) 

Since, by (35), /3(ra) - 0(m - x) < \2 a{m \ then (36) holds for 2 a °" ) ^ 
m < 3[2 ]. Substituting these values of a(x), 0(x), a(m — x) and 

/3(ra — x) into (32) and using (24) gives 

Fl . (m) _ a(m) (^£-) + 4 — >(L=^), (37) 

which completes the induction for Case 1. Conversely, if 2 ai '" ) ^ m < 3 
[2 a m)_L ] then .ro = 2° '" _1 and, for any x with .r ^ x ^ m — .r , we have 
a(x) = a(ra — x). Hence (32) is minimized for x ma x = xo = 2 a( '" ) ~ i by 
the above argument. 

Case 2. Since we are now considering only possibilities outside Case 1, 
we must have 3[2 a(m)_1 ] g m ^ 2 a(m)+1 . Using (29), and the fact that 
Xo ^ x ^ m — x , it follows that either a(x) = a(m — x — 1) = 
a(m — x) — 1 , so that m — x is a power of 2 and x < ra/2, or else 
a(m — x) = a(x — 1) = a(x) — 1, so that x is a power of 2 and x > ra/2. 
In the former alternative, we obtain 

ra = x + (m - x) = 3[2 oU> ] + p(x), 
a(m) = a(x) + 1 and |8(w) = 2 a(x) + /3(.r) = x. It follows that 

•r = x a „ = m - 2 a(m) . (38) 
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Iii the latter alternative, it can be concluded that x = 2" T,i > m/2. In 
the latter alternative (or in both alternatives), we now compare the 
values of Ci*(x), computed from (24), for the two arguments, 

Xl = 2 a{m) > m/2 and x 2 = m - 2 a(ra) < m/2. 

Using the general result that p(m*) < m*/2 for any positive integer m* 
and the implied result that a(m - 2 a(m) ) = a(m) - 1, it is easy to show 
(the details are omitted) that Ci*0r 2 ) < &*(*i) for < q < 1. Hence, 
for either alternative, we obtain the same result, (38). As a corollary, 
(22) holds. In this case Ci*(x), after algebraic simplification, is given by 

CV(s) = [aim) - 1] f^Q + r*» [4^] » (39) 

which, using (24), again gives the result (23). The fact that (23) holds 
for m = 2, 3 and 4 is easily shown and the details are omitted. This 
completes the proof of (23). 

It is also possible to get expressions for G t *(m, m), Gi*(m, n) and 
Hi(n) for large values of q using (23) and the fact that, for any posi- 
tive n and sufficiently large values of q (depending on n), 

H x (n) = 1 4- (1 - q n )Gi(n, ») - 1 + p&*(», n). (40) 

Then, by (18), we obtain for large values of q 

G!*(m, m) = m - 1 + qFfim) 4- p £j Fi*(j) t (41) 



J = 2 



Gi*(m, n) =n-l-(n-m- l)q m 4- F x *(m) 

n-m-l n-1 (42) 

+ p(l " q m ) E Fi*U) + V E Fx*{j), 



j=2 j=n—m 

n-1 



fft(n) = q + np + pF!*(n) + p 2 E ^*(i), (43) 



j-2 



where (42) is to be used only for m < n and the summation from a to b 
is taken to be zero for b < a. Using these and (23), the last equation 
(i.e., the equation for the last g-interval which ends at unity) can be 
obtained independently and more simply for Fi*(m), Gi*(m, n) and 

Hi(»). 

It is clear from the results above that, as q approaches unity, 

limffi(n) = 1 (43a) 

and that, for q in an interval ending at unity, we have x H (n; q) = n. 
Here x„{n; q) is the size of the next test group when the procedure R^ 
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is used for an H(n) -situation with a priori probability q of a unit being 
good. In an //-situation the probability that there are no defectives pres- 
ent approaches unity as q approaches unity. In a G-situation with m > 1 , 
the probability that there is exactly one defective in the defective set 
and none in the binomial set approaches unity as q approaches unity. 
Defining x g (m\ q) in a similar manner, it was shown above that 
x a (m; q) = x,„ nx [as given in (21)] in an interval ending at unity, and 
it follows from (23), (41) and (42) [or (23) and (15)] that, as q ap- 
proaches unity, 

lim Fdm) = a(m) + 2 ^ (m > 1), (43b) 



lim G\(m, n) = 



I + lim Fi(m) = 1 + a(rn) + 2 ^^ (for n > m) 



m 



(43c) 



1 - - 1 + a(m) + 2 ^ (for n = m). 



m m 

It is interesting to note that the result (43c) depends on m but not on n. 

VII. CONJECTURED PROPERTIES OP Ri 

In this section we shall state some properties which appear to hold 
for procedure 7?i based on numerical calculations for N ^ 16 but have 
not been proved for all N. 

1. For any G-situation with fixed m > 1, if rr G (m; q) denotes the size 
of the next test group under R h then x aim; q) is a nondecreasing step 
function of q with step size unity. That is, for any pair q ^ q + 2£ 1, 

Xaim; q) ^ :r (ra; q + ) (44) 

and, for sufficiently small e > 0, 

xaim; q + e) ^ x a (m; q) 4- 1. (45) 

Also for fixed q the value of x a (m + 1 ; q) is either the same or one greater 
than Xa{m\ q); i. e. 

xaim;q) ^ x a (m 4- 1; q) ^ x a (m; q) + 1. (46) 

If the dividing point between x and x + 1 for any (/-situation under R x 
[denoted by q Q ix\ m)} is shown to exist (and be unique), then the three 
properties (44), (45) and (46) are equivalent to the two properties that, 
for any m > 1, 
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q a (x; m) ^ q (x; m + 1) for 1 ^ x < x mBX , (47) 

q G (x; m) < q (x 4- 1; m) for 1 ^ x < .r innx - 1. (48) 

The assumption used in Section VI that the largest rvalues are as- 
sociated with the largest q- values is a simple consequence of (44). 

2. For any //-situtation, we can define x H (n; q) similar to x (m; q) 
and the property corresponding to (44) still holds: that, for q ^ q + ^ 1 
and all positive integers n, 

x fl (n;q) ^ x H (n;q + ). (49) 

We can also define q H (x; n) similar to q a (x; m) with the understanding 
that, if x does not appear under Hi{n) in Fig. 3, then the interval for 
x is assumed to have length zero but the endpoints still exist, and, in 
fact, q H (x - 1, n) will then be equal to q H (x; n). Then (49) is equiva- 
lent to the property that for all x ^ 1 and all positive integers n 

qn{x;n) ^ q H (x + l;n). (50) 

3. If the experimental situation is such that it is impossible or eco- 
nomically impractical to identify or keep separate the individual units 
in any test group, then, after each test on a batch of x units, the disposi- 
tion of the x units must be made on a batch basis. In such a situation 
it is conjectured that the procedure Ri is the optimal procedure for all 
values of q. 

4. There are various patterns existing in Table II B both within a 
column and across columns, none of which have been proved. For ex- 
ample, if q = q a (x; m) is the dividing point between x and x + 1, then 
the first entry in the appropriate column is 1 — q 2x — q x \ and the last 
entry, f or m = <», is 1 — q x — q I+1 . In the first column, the entry can 
be written as 1 — q — q + q" for all m. In the second column of Table 
II B the pattern shown is to replace the highest power of q (say, q) by 
q h ~ 2 or by q h+2 + q h+3 , depending on whether m is odd or even. Thus the 
pattern displays a cycle of 2. In the third column one can similarly find 
a pattern with a cycle of 3, starting with m ^ 9. If a general rule for all 
these patterns were proved then it might be easier to find the dividing 
points for higher values of m. The conjecture in this case lies in the fact 
that these patterns exist and can be mathematically established. 

VIII. CHARACTER OF Ri FOR SMALL VALUES OF q 

For the procedure Ri it will now be shown that, when 
q < qo = KV5 - 1) = 0.G18 
(to three decimal places), then, for both G- and //-situations with any 
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positive integers m ^ n, the units are all tested one at a time. Of course, 
if we start with the //-situation and test units one at a time, then a (/-situ- 
ation never arises, but in the induction proof that follows it must first be 
shown for the G-situation and then for the //-situation. This property 
that units are tested one at a time was recently shown 4 to hold for the 
optimal procedure (without specifying what the optimal procedure is or 
whether it exists). Simple formulae for H x (n), Gi(m, n) and F\{m) are 
obtained for q ^ q . 

Theorem 2: For procedure i?i with 1 :g m ^ n and ^ q < q , 

x a {m; q) = x H (n; q) = 1, (51) 

Hi(n) = n, (52) 

m-1 

Gi(m, n) = n - ™_ m , (53) 



»— 1 



ri (m)-i+ l -*~ Z mm . (54) 

V 1-5"* 

[Remarks: The last term in (53) results from the possibility of saving 
one test if the defective set of size m contains exactly one defective unit 
that is discovered inferentially by showing that all the other m — 1 units 
are good. It is interesting to note that (54) can be obtained by sum- 
ming the series 

F »( OT ) = r^i (1 + 2ry + 3r/ + • • • + (m - 1)^ 

1 " q (55) 

4- (m - Dq"- 1 ). 

In the proof below, (54) is shown first and then (52) ; the proof of these 
contain the result (51). Then (53) follows from (15), (52) and (54)]. 
Proof: The proof of (54) is by induction. The result holds for m = 1, 
since /'(l) = 0. Assuming (54) holds for arguments less than m, we 
can use (19) with (16) to obtain 

Fl *( m) = LzX + ^ lAl V" 



• /-I -1 

mm \q I - 



V J 

+ I - q*~ l - (x - l)pq x \ m p - g*"- 1 - ( m - Ijpf l (5fi) 

+ ~ir- + p ,jSl |x(4 '" ~ ' n ~ m "' ] - 
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To prove part of (51) it is now shown that the minimum of the expres- 
sion above in braces [say, f (x)] is attained at x = 1 . Since this is obvious 
f or m = 2, it is now assumed that m ^ 3. Then 

fo(x) - Ml) = (x - \)q m + 1 - xq* - pq*~\ ( 57 ) 

and it suffices to show that, for q ^ q and X = 2, 3, • • • , m — 1, 

f x (x) = 1 - xq x - pq 1 ' 1 ^ 0. (58) 

Similarly, it suffices to show that, for q ^ go and x = 3, 4, • ■ ■ , m — 1, 

ff x ) = *<*> ~ - fl(2) = 1 + q - xq*- 1 - pq*- 2 ^ 0. (59) 
Q 

More generally, if it suffices to show that, for q g q and x = y + 1, 
y + 2, ••• ,m - 1, 

fv(x) = 1 + (y _ l)g - .rr/-" +1 - pr/"" ^ 0, (60) 

then it also suffices to show that, for q ^ q and a; — y + 2, y + 3, • • ■ , 

m — 1, 

w(x) = A(g) -/»>(?/ + D = ! + yg _ ^ _ p^r-i ;> . (61 ) 

Setting?/ = m — 3 in (61), it suffices to show that, for q ^ q and x = 
m — 1 (we can now replace x by m — 1), 

1 + (m - A)q - (m - 2)q* = 1 - q - q + (m - 3)pg ^ 0. (62) 

Since q is the root of 1 - q - q and m ^ 3, the inequality (62) is proved 
and the minimum of (56) is attained at x = 1. Setting a; = 1 in (56) 
gives the bracketed expression in (56) and proves the result (54). 
To prove (52), we substitute (54) in (20) to obtain 



(63) 



Hi{n) = 1 + min \q x (n - x) - pq* ' + pE H % l 

I p I (n - Off* -1 
,-i J 

= n + 1 - max {xq' + pg 1-1 }, 



l<x<n 



where the value of ?* ! for q = and .r = 1 is taken to be unity. 

To prove the rest of (51), it is now shown that the maximum of the 
expression above in braces [say, h(x)] is attained at x = 1 for q ^ q . 
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Then, for q ^ qo and x ^ 1 , 

h(x) — h(x + 1) = xpq x — 2q z + q x ~ l 

(64) 

= q [l - ff - e + (« - l)pg] ^ o. 

Clearly, the maximum of /i(.r) is attained at x = 1 not only for 
< 9 ^ #o but also for g = 0. Setting x = 1 in (63) gives Hy(n) = n, 
and this completes the proof of (52). The fact that the minimum is 
attained only at.r = 1 for q < q in both (56) and (63) proves (51) and 
shows that, under Ri with q < q Q , units are tested one at a time. 

IX. CONSTRUCTION OF TABLES FOR Ri 

Figs. 3 and 4 describe the procedure Ri for n = 2(1)16 and m ^ n 
in the form of two diagrams that are easy to use in a practical situation. 
Tables III A and III B give the polynomials, the roots of which are the 
dividing points in Figs. 3 and 4. Table IV A gives the polynomial equa- 
tion for Fi*(m) for m = 2(1)16. Tables IV B and IV C give the poly- 
nomial equations for Hi(n) and G*(m, n), respectively, for n = 2(1)12 
and 2 ^ m ^ n. These can be obtained from (6), (7) and (8), or from 
(18), (19) and (20) and the boundary conditions, H x (0) = Fi*(l) = 0. 
For the sake of brevity, Table IV C has been reduced so that it gives the 
results only for q ^ 0.85, and only for pairs m,n which can arise starting 
from an H(n) -situation with n ^ 12. 

Having computed Hi(n) and F*{m) for 2 ^ m :§j n ^ 12, we can 
make the procedure Ri explicit for 12 ^ n ^ 16 by a different method, 
which will now be explained. Let H}(n \ x) denote the value of Hi{n) if 
(i.e., for those q- values for which) the next sample size is re; let F*(m \ x) 
be defined similarly. Then (20) can be written as 

tfi(n | x) - Hxin - 1) = 1 + pFSix) 

^ , (65) 

- E q[Hi(n - i) - Hrin - i - 1)]. 
t=i 

Writing a similar equation for H\{n \ y) for y > x and subtracting gives 
ffi(n | x) - H x (n | y) = -p[F,*(</) - F,*{x)] 

E q\H,(n-j) - H^n-j - 1)]. 



(60) 



1=1 



If 4 ^ x < y ^ 16 and 12 < n ^ 16, then the right member of (66) 
involves //-function arguments only up to 12. In particular, for y = 
x + 1 we set the right member of (66) equal to zero and obtain a poly- 
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nomial whose root (between zero and one) is the dividing point, q H (x, n), 
between x and x + 1. Table II A shows that, for n > 12 and x < 4, the 
pattern for the dividing points is well stabilized; for example, the value 
of q H (l, n) = 0.618 (to three decimal places) is shown to hold for all 
n in Section VIII. By considering various pairs x, y (most often of the 
type x, x + 1) in (66) it is possible to determine the procedure Ri for 
the H(n) -situation for 12 < n ^ 16, without explicitly computing the 
formula for H^n). 

Similarly, we can do the same for the G-situation by using 

Fi*(m | x) - F 1 *(m \ y) = F 1 *(y) - F x *(x) 

( 67 ) 
+ tfFfim - y) - q'Ffim - x), 

but in this case (according to Section VI) we need only consider values 
y = x + 1 up to and including m/2. 

Table II A gives a numerical comparison for N = 4, 8 and 12 of R\ 
and several other procedures, two of which are based on the work of 
Dorfman 1 and Sterrett; 3 the others are defined in the Appendices to this 
paper. Table II B gives a brief numerical comparison of Hi(n) and H 3 (n) 
(corresponding to procedure R 3 denned in Appendix B) for large values 
of n [viz., n = 10(10)100] and q = 0.90, 0.95 and 0.99; these entries 
were computed on the IBM-704. 

Tables V A, V B and V C give the numerical values of H x {n) and 
Gi(m, n) as well as the values of x {m; q) and x H (n; q) for q = 0.90, 0.95 
0.99 for 2 ^ n ^ 100, and for appropriate values of m; these entries were 
computed on the IBM-704. For q = 0.90 the value of x is always at most 
9 and hence, if we start with an //-situation, there is no need to consider 
values of m > 9; similarly, for q = 0.95 we disregard values of m > 19. 
For q = 0.99 we should consider all values of m up to and including 
m = 100 but many of these were omitted for the sake of brevity. It is 
interesting to note that Gi(m, n) is strictly monotonic in the second argu- 
ment (and hence also in the argument n — m) for fixed m, but it is curious 
and difficult to explain why it is not monotonic in the first argument for 
fixed n. 

X. A SUGGESTED PROCEDURE FOR THE CASE OF UNKNOWN q 

It is reasonable to expect that a knowledge of good procedures for the 
case of known q will suggest good procedures for the case of unknown q. 
From this point of view we consider modifications of the basic procedure 
Ri that make it adaptable when q is unknown. It is suggested that after 
each test we form a new estimate of q and that the procedure R\ be 



GROUP TESTING TO ELIMINATE DEFECTIVES 1201 

used with the estimated value in place of the true value. At the outset 
we can start with an estimate based on past experience or we can start 
by testing one unit at a time. A thorough investigation of the relative 
merit of this procedure has not been carried out. Some discussion on the 
maximum likelihood method of estimating q is given below. 

Let d and s denote the number of units proven defective and proven 
good, respectively, so that at any stage of experimentation we have 

N = d + s + m + (n - m) = d + s + n, (68) 

where N is the total number of units at the outset, m is the size of the 
defective set (which is known to contain at least one defective) and 
n — m is the size of the binomial set. The likelihood L of the observed 
result (68) is given by 



(N ~n\ d 



. N-n-d 



L=" H (1 " <f). (69) 



Then it is easily shown that 

f-(logL) = -± (logL) = -- 

i/t/ dp pq _ 



d- (N -n)p+,pW- \. (70) 



mpq J 



Setting the latter equal to zero, we find that, for m ^ 0, the maximum 
likelihood estimate q of q is the root of 

(JV - n + m)(l - q")(l - q) - d{\ - q") - m(l - q) = (71) 

or, equivalently, 

m 

s-dE (/ - (m + s)q" = (72) 

and, for m = 0, we have q = s/(d + s), the usual estimate. For s = 
and m + d ^ 1, we get q = and, for s = 1, it is easily seen, using the 
Descartes Rule of Signs, that (72) has exactly one root q (allowing multi- 
plicities) in the unit interval and hence q is uniquely defined. The re- 
maining case, s = m = d = 0, can only occur at the outset when there 
is no observations on which to base an estimate. It is interesting to note 
that the same result (71) or (72) can also be obtained by computing 
the conditional expected proportion of defectives among the N units, 
given the observed s, d, m and n, and setting it equal to 1 — q. The 
equation thus obtained is the same as (71), and its root is q. 

The above method of getting an estimate is being suggested in con- 
nection with procedure Ri , but it can also be used in connection with 
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the procedures R 2 and R A (see Appendices) without any change. For 
procedures #3 and R 6 we can have several defective sets and several bi- 
nomial sets at any one time, and (71) then becomes 



d+(l-q) 



±T^f--( N -p'')] = °- (73) 



where mi , m 2 , • • • , mi are the sizes of the defective sets, n/, w 2 ', • ■ ■ > n -> 
are the sizes of the binomial sets, and 

/ j 

N = d + «+E»*+2>/- ( 74 ) 

If the number of tests carried out is large (and hence N — n must be 
large), the maximum likelihood estimate is approximately normally dis- 
tributed with expectation equal to the true value of q and variance 
given by 



*»« [*(*$*)!• (75) 

where m and n in (69) are to be regarded as chance variables. Taking 
expectation first for fixed m and n and then with respect to m and n, 
gives 

*»"[ J (T) + *(r£)T- (70) 

Since mpq m ~ 1 < 1 — q m for all wi and all q < 1, and since the expecta- 
tion of a square is nonnegative, for asymptotically large N — n 



2 



g_g < */an < Pg (77) 

p + gF - E(n)\ ~ Kq ' ~ N - E(n) ' K 

In particular, if we continue to test until a fixed proportion > of 
the N units are determined to be good or bad (i.e., until n/N =1 — 6, 
approximately) then, for asymptotically large N (so that N6 is also 
large), we obtain 

P 2 9 <^( r ;)<gg (78) 

Ndp + q- W/ ~ N6 

For 6 = 1 and large N, the two bounds are essentially equal and the 
common value is the same as for ordinary binomial sampling. In general, 
at any stage of experimentation it appears to be conservative to esti- 
mate a 2 (q) in the same way as for ordinary binomial sampling based on 
N — n observations, using the value of n that is actually realized at that 
time. 
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In regard to the procedure, if the size of the very first test group is 
based on past experience, the question arises as to whether this past 
experience should also enter into the second, third and other early esti- 
mates of q. If it does not enter, then in the early tests we may find sud- 
den jumps from testing very small numbers to testing very large num- 
bers and vice versa, both of which are undesirable. This makes it useful 
to find a method to continue to use past experience until the estimate of q 
(without using past experience) is stabilized. In the absence of past ex- 
perience, this same feature may make it desirable to test several units 
one at a time before starting to use any group-testing procedure. 

XI. SOME GENERALIZATIONS OP Ri 

Returning to the case of known probabilities q, we consider some gen- 
eralizations of the same basic problem and, in each case, the appropriate 
generalization of the procedure Ri . The appropriate formulae will be 
given, but only a few simple computations will be carried out. 

1. Two (or more) different kinds of units with known probabilities 
(say, qi ^ q 2 ) of a good unit are present and both can be put into the 
same test group. 

Let H u (ni , n 2 ) denote the expected number of tests required under 
the proposed procedure Ru if there are n< units of type i with q = 
qi{i = 1,2) and the binomial chance variables associated with the units 
arc mutually independent. Let G u {m 1 ,m 2 ;ni, n 2 ) denote the expected 
number of tests required under R n if there is a defective set containing 
mi units of type 1 and m 2 units of type 2 (known to contain at least one 
defective among the mi + m 2 units) and a binomial set containing 
m — mi ^ of type 1 and n 2 — m 2 ^ of type 2. The recursion formulae 
corresponding to (G), (7) and (8) are 

Hn{n\ , n 2 ) = 1 + min { qi x q 2 v H n («i - x, n 2 - y) 

(7y) 
+ (1 - qi x q 2 v )G n (x, y;ni,n 2 )\, 

where the minimum is over pairs (x, y) with ^ x ^ Wi , ^ y ^ n 2 
and x + V ^ 1 ? and 

{/ x V "'l/t m 2\ 

("hr ml m! ) 
\ l - qi <?2 ' / 

•Gn(wii— x, m 2 — y;n x — x, n 2 — y) (80) 

, / i - gi!gg^ \ .1 

+ ( 1 _ »., »- 2 )Gii(x, y; m , n 2 ) > , 
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where the minimum is over pairs (x, y) with ^ x ^ mi , ^ y ^ ra 2 
and 1 ^ x + 7/ ^ mi + m 2 — 1. The boundary conditions state that, 
for all (ft £i <fe > 

Gu(l, 0; n, , n 2 ) = ff u (m - 1, na) for all n t £ 1, n 2 > 0, (81) 

(?n(0, 1; ni , n 2 ) = ff u («i , n 2 - 1) for all n, ^ 0, rc 2 ^ 1, (82) 

Gn(mi , 0; «i , 0) = Gi(m, , n, ; gi) ; 

Gn(0, m 2 ; 0, n 2 ) = Gi(m 2 , n 2 ; q 2 ), 

ffn(ni , 0) = T/^n, ; ft ); ff„(0, n 2 ) = H^n, ; ry 2 ), (84) 



(83) 



where the right-hand member of each equality in (83) and (84) refers 
to the basic procedure Ri denned by (G), (7) and (8). 

It is clear that tf u (l, 0) = ff u (0, 1) = 1 and H n (0, 0) = 0. It fol- 
lows from (80) that, for qi ^ q 2 , 

G 11 (l,l;l,D = 2 T g2 ~ glg2 > (85) 

1 - tfl?2 

and the rule is to test first the unit of type 2. Using this result, we can 
compute 

H n (h 1) = 1 + min(l, 2 - q 2 - Ql q 2 ), (86) 

and the rule is to test either unit separately if 1 — q 2 — q\q% > and to 
test both simultaneously if 1 — q 2 — q\q 2 S 0. The latter inequality is a 
direct generalization of the inequality 1 — q — q 1 ^ 0, which played a 
prominent role in the basic procedure R\ . We state (without proof) that, 
if q\ ^ #2 < \ (\/5 — l) = 0.G18 (to three decimals), all testing is car- 
ried out one unit at a time. 

2. Two (or more) experimenters may be working on a single set of N 
units by carrying out simultaneous, parallel group-tests and cooperating 
in such a way as to minimize the time required to accomplish the task. 

It is clear that no saving can be effected in the expected total number 
of tests by having more than one experimenter. However, if the simul- 
taneous tests are regarded as a stage, each of which lasts the same amount 
of time, then minimizing the expected number of stages is equivalent 
to minimizing the expected time required to accomplish the task. These 
remarks indicate that there may be some conflict in these two aims of 
reducing the expected time and the expected total number of tests. For 
this reason, it should be stated that our primary emphasis in this prob- 
lem is to reduce the expected time. 
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Let m and m denote the sizes of defective sets and let n — ( w + in) ^ 
denote the size of the binomial set. Let H n (n) denote the expected num- 
ber of stages required for m = m' = by the proposed procedure R i2 . 
Let G v <(m, m', n) denote the expected number of stages required by R vl 
if we have two defective sets of size m, m' and one binomial set of size 
n — m — m' ^ 0. Let Gi 2 (m, 0, n) and Gi 2 (0, m, n) be denoted by 
G i2 (m, n), so that G, 2 (0, n) = H i2 (n). The recursion formulae for R i2 are 

ffii(n) = 1 + min [q x+u H 12 (n - x - y) + f/(l - q v )Gi 2 (y,n - x) 

■r. a ^ 1 

s + U^n (87) 

+ r/(l - q*)Gu(x,n - y) + (1 - r/)(l - q") G 12 (x,y,n)}, 

1 g u g "»'-i 

■Gn(m — x,?n — y,n — x — y) 

+ (f^")(^)^ -**.-»> (88) 



+ 



and 



:s<j) 



Gn(m,n) = 1 + min / min Gi 2 '(a:, //), min G i2 " (x, y) 

where Gu(x, y) and Gn {x, y) are defined by 

Gi 2 '(.r, </) = \ x _ m ) Gm(»i - x - y, n - x - y) 

1 5 (90) 
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and 

Gli "( Xf y ) = g Y_~; ; Gnim -x,n- x - y) 

+ g!<l ~ p Gl2 (x, n - y) 

1 - « (91) 

(l-gW-<T) (? 12 ( W _ *, ,, n - x) 
1 - q m 

1 — q m 



(92) 



The boundary conditions state that, for all q, 
G l2 (m, 1, n) = G , i»(l, m, n) = G l2 (m, n - 1) 

for ^ m ^ n — 1, 

Gu(l,n) = ff«(n - 1) for n ^ 1, (93) 

tf 12 (0) = 0. (94) 

It is easy to see that H l2 (\) = G u (2, 2) = H l2 (2) = 1 and 
G, 2 (3,4) = G 12 (4, 4) = 2 for all g. 

Remark 6: The extra complication in (89) insures that, for n ^ 2, 
one experimenter will not be idle while another is carrying out a test. 

Remark 7: It is conjectured that in (89) the possibility x + y = m 
can be omitted, with the exception of the single case m = n = 2 (and 
m' = 0). 

Remark 8: It is also conjectured that G 12 ", which is needed for the 
cases n > m = 2, can be disregarded when m > 2; i.e., that G i2 always 
gives a smaller minimum for m > 2. 

Remark 9: It is conjectured that, at any stage in which m = m' ^ 2 
or in which we have both m = m' = and n even, the two test group 
sizes, x and y, will be equal. If either m or m' = 0, it is conjectured that 
the two test group sizes will differ by at most unity. 

Further calculations yield 

Test Group Sizes 
x y 

Gl2 ( 2) 3) - \±1 for a11 ? 1D 1B > (95) 

1 + 9 

Gf 12 ( 3j 3 ) = 2 + 2< * + g' for all q ID ID, (96) 

1 + q + q 2 
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(2 forO ^ q < 0.682 ID 15, 



Gi 2 (2, 4) = 



3 + q ~ q 



1 +5 



(97) 
for 0.682 ^ q ^ 1.000 ID 25, 



G 12 (2, 2, 4) = 2 ± g +/ for all q ID ID', (98) 

(i + q)- 

(2 for ^ g < 0.691 15 15, 

ffn(4) = (99) 

1 3 - 3g 2 + 2g 3 - q 4 for 0.691 ^ q ^ 1.000 25 25, 



where 15 indicates that 1 unit is taken from the binomial set to form 
one of the two test-groups and D, D' denote different defective sets. 

It is interesting to compare the above result for i/i 2 (4) for q ^ 0.691 
with the procedure R\* of giving each experimenter two units to analyze 
independently of each other and without any mutual cooperation. Let T 
denote the total number of tests and S denote the number of stages re- 
quired. Then, for q ^ 0.691 (letting Ti , T 2 denote the number of tests 
in two independent experiments with n = 2 under 5i), it is easily shown 
that 

E{S\ R n *\ = £{max(T, , T 2 ) | RJ = 3 - q 2 - q 4 , (100) 

E\T | R n ] = 2// 12 (4) - 2q - 2q~ 4- 2q - 2q\ (101) 

Hence we find that, for q ^ 0.691, 

EiSlRn) - E{S\R n *} 

= // 12 (4) - (3 - f/ - q) = -2r/(l - q) g 0, 

E\T\R n ) - E\T\R n *\ 

= 6 - 2q - 2(f + 2q 3 - 2q - 2H l {2) = 2q\\ - q) ^ 0, 



(102) 



(103) 



which illustrates the fact that /? J2 effects an improvement in the expected 
number of stages at the expense of a slight increase in the expected total 
number of tests. 

3. In this generalization we apply the restriction that any one unit 
can be included in at most K group-tests. This is particularly appropriate 
in the blood testing application, where a single blood sample can be used 
in a small number K of blood tests and the patient does not want to be 
annoyed by having more than one blood sample taken. 

In this problem there is again only one defective set but it is now de- 
noted by a vector m = \m , ??h , • • • , itik-i} , where m } - ^ is the nura- 
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ber of units that have already been included in j group-tests. Similarly, 
the union of the binomial and defective sets is denoted by 

n = {n , fti , • • • , riK-t} 

and the binomial set alone is the difference n — m. The symbol for 
the size of the next test group will be x = {x , x% , • • • , x K -i\, where Xj 
is the number of units taken from n, in an //-situation (from m,j in a 
G-situation). The symbols x, m, n will be used for the sum of the com- 
ponents in the vectors x, m, n, respectively. Let Gn (m; n) denote the 
expected number of group-tests required to remove all defective units 
if the defective set is m and the binomial set is n — m. If m = 0, 
we denote this expectation by H 13 K (n). For the special case in which n 
has all except one component (say, rij) equal to zero, we will drop the 
zeros and write H i3 K {nj), with a scalar argument. The recursion formulae 
for this procedure R U K are given by 

H 13 K (m) = 1 + min {q x H 13 K (n - x) + (1 - q*)G u K (x; n)} , 

all x with (104) 

lgi^n and N ' 

g x, < nj 
W-0,1, .. .,*-!) 



Gu'dn; n) - 1 + min \(t £) G 13 K (m - x; n - x) 

all x with ^V — Q / 



1 gigm-l 

< x: < mi 

(/ = (), ••■■K-D 



(105) 



i[^]^x ;il i 



where, as usual, m > 1. The boundary conditions state that, for all q, 
we have Hu K (0) = and, for m = 1, we can write Gi 3 A (m; n) = 
Hu K (n — m). It is easy to see that His (n) — 1 for n — 1 and all q. 
Some further computations for n = 2 and n = 3 give, for any K: 

q -interval a; -value j- value 

(ffi(2) ^ q g 1.000 (see #i) ^ / ^ K - 2 (106) 

\2 g q ^ 1.000 (0, • • -, 0, 0, 1) j = K - I 

ffl3*(3y) = 

^ q g 1.000 (see Ri) ^ j £ K - 3 

Org g =£0.618 (0, •• -,0,1,0)1 (107) 

0.618 g q ^ 0.707 (0, • • •, 0, 1, 0) \ j = K - 2 

0.707 g?^ 1.000 (0, •••,0,3,0)J 

^ g g 1.000 (0, •• -,0,0, 1) j = K - 1. 
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If K = 1, all units are tested individually. If K = 2, then, after a set 
of size m > 1 is shown to be defective, all the units in that set are tested 
individually; this is the procedure recommended by Dorfman. 1 A more 
thorough investigation for particular values of K ^ 3 will be given in a 
separate paper. 

Under R\ for any //(n) -situation the maximum number of additional 
tests M*(n) in which any particular unit will be included before experi- 
mentation is concluded (allowing a random reordering of units in the 
binomial and defective sets after every test) occurs when q is close to 
unity and the sample has all units defective. It is easily seen that, un- 
der R\ , 

M*(n) = M(n) - (to - 1) = (to + l)a(n) + 3 - 2 1+o( " ) , (108) 

where a(n) and M(n) are defined by (10) and (11). Under R ri K , let 
,/(n) denote the largest subscript associated with a nonzero component 
of n. For K ^ M*(n) + j(n), the restriction that no unit should be in- 
cluded in more than K group-tests does not affect the procedure Ri , 
and hence we have for procedure Rn K 

H ls K (n) = Hi(n) for ^ j(n) ^ K - M*(n), (109) 

which generalizes some results in (106) and (107) and shows that R ri K 
is a generalization of R t . 

Since units in the last component of the defective or the binomial sets 
cannot be tested in groups, we can remove at any time for individual 
testing all units in the last component of the binomial set and all but 
one of the units in the last component of the defective set without af- 
fecting the expected number of tests under R l3 K . It is easy to show that 
this leads to the two reduction formulae 

// 13 A '(n) = riK-i + //./(n , • • • , to*_ 2 , 0), (110) 

(\ - n -1+m *- 1 \ 
G 13 *(m; n) = n K -i + ( — x _ m ) ^« («o . • ' ' . ^- 2 , 0) 

• [Gn K (mo , • • • , m K -2 , 1 ; n , • • • n K -z ,0) — 1], 

which are useful in computations and for checking. 

It is conjectured that, for m = or m > 1, the procedure Ri 3 K can 
always be carried out by putting in the next test-group only units that 
have been included in the same number of group tests (i.e., units in the 
same subset) ; the only possible exception to this is that, in any //-situa- 
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tion, if q is sufficiently close to unity, then R U K will call for a test of all 
the remaining units; i.e., x = n. 

Under the above conjecture, it is possible to carry out a simplification 
as in Section V and show, in direct analogy with (19), that 

m 

Fu K (m) =Ef' + min lq x Fn K (m - x) + F 13 **(x)), (112) 

where x has x { in the (i + l)th position and zeros elsewhere (so that 
x = Xi), and i is defined as the subscript associated with the first non- 
zero component of m. The function Fi 3 *(m) is defined as the expected 
number of group-tests required to reach the next //-situation, and then 
we define, as in Section V, 

Fu K (m) =(L=-£ S )F 13 K (m). (113) 

Hence under the above conjecture it is again seen that, for any G-situa- 
tion with m > 1, the next test-group x depends on m but is independent 
of n — m. 

It appears to be true (but has not been rigorously proved) that, in 
this case also, for q < q = 0.618 (to three decimal places) all units are 
tested one at a time. 

XII. AN ASYMPTOTIC FORMULA FOR H\{n) 

In this section we shall use results obtained by considering an informa- 
tion procedure R 2 , which is defined in Appendix A. The procedure R 2 
appears to be a best test in the sense that it maximizes the information 
in the very next test but does not take into account the exact finite num- 
ber of units present and the possible ways of distributing them among 
subsequent tests. It is therefore intuitively reasonable to expect that the 
procedure it!i tends toward R 2 in the ^-situation as n -> » and also in 
the (r-situation asm-* °o . A more rigorous proof of this assertion would 
be desirable. It should also be pointed out that there is considerable 
numerical evidence in Tables IIIA and IIIB of the above assertion, which 
explains the reason for putting opposite m = °° and n = °° in these 
tables the polynomial equations 

1 _ q x -q x+l = {x= 1,2, ••• ), (114) 

which are derived for procedure R 2 in Appendix A. 

We shall now derive an asymptotic formula for Hi(n) for large n 
based on the assumption that the above reasoning is correct. For large 
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values of n and fixed q, the expected number of tests required under 
procedure R x is approximately given by 

(expected number of tests needed to reach \ 
the next //-situation under Ri\ , . . _ .. 
: . (115) 
expected number of units analyzed between / 
H -situations under Rj 

The ratio of n to the denominator in (115) is the approximate number of 
.//-situations reached if we start with n units, and this is clearly to be 
multiplied by the expected number of tests required to proceed from one 
//-situation to the next. Let T and U denote the chance variables in the 
numerator and demoninator, respectively, of (115). For a fixed q we 
find from the limiting procedure R2 that, for an H(n) -situation with n 
large, we will, under Ri , "almost always" be using the same test-group 
size x, where x is that positive integer for which q x is closer to one-half 
than is either q x ~ l or q x . Then, for this fixed integer x, which depends 
on the given q, we have 

E{U I RJ = xq 1 + V ZjV'- 1 = \^i- , (116) 

;=i 1 — q 

which is obtained by assuming a single randomization of the order of the 
units at the outset and considering the different possible positions of the 
first defective. 

Since Fi(x) is the expected number of tests required under Ri to get 
from a G(m, w)-situation to the next H{n) -situation, we have 

E[T I fij = q* + (1 - r/)[l + F^x)] = 1 + V F 1 *{x), (117) 

where F*{x) is tabulated in Table IV A for x — 2(1)16 and all values of 
q. Hence, we obtain from (115), (116) and (117) 

HM S- npil ± VF1 * (X)] , (118) 

1 - q* 

where x is defined above in terms of q. This is the main result of this 
section; we now consider some special cases. 

For values of q close to unity we can use (23) to replace F*{x) by an 
explicit expression. If we also replace (1 — q*)/{\ — q) by z for q close 
to unity, Ave obtain for q close to unity 

HM «- {1 + p[xa(x) + 20(a)]}, (119) 

x 

where a(x) and 0(x) are defined in (20a). If q approaches unity, x be- 
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comes large; if a(x) = 2, then 20 (3) < 2 a(x)+x ^ xa(x) and, since 
a(x) — > », it follows that we can disregard 2/3(x) in (119). For q close 
to unity and x large, the dividing points get closer and closer, and we 
obtain 



so that x = 



so that 



log 2 1 - 



S 1 - <1 ~ q = 1 - 2f/, 
-1 

. Also, from the definition of a(x), 



Q/J 



2 aw < x < 2-2 



a(x) 



log 2 ( ^ J = a(.r) = log 2 (x). 



(120) 

(121) 
(122) 



Using the upper value in (122) gives 



Q + np loft [log s Q] . (123) 



Ih(n) ^ n log 2 ( - ) + np log 2 

It can be shown that the first term in (123) goes to zero faster than the 
second as q approaches unity, and hence we drop the first term and re- 
write the second in the form 

H x (n) S -nplog 2 p. (124) 

In particular, if p = 1/w so that 2=1— (l/n), we obtain from (124) 

Hi(n) ^log 2 n. (125) 

In this case, we also have x = n — 1/p, and a better estimate is obtained 
by setting q = e~ l in (23) and (118). Assuming that (23) is either equal 
to Fi*(x) or is a good approximation to it, we obtain 



Hi(n) ^ - — zr — = - 

x _ 1 — e 1 J x 



1 + a(x) + 



,2/3(x)/x-| 



e - 1 



. (126) 



For the case n = N = 100 and q = 0.99, we obtain 8.32 as the exact 
value from Table VC; log 2 100 = 6.64, from (125); 7.79, from (119); 
and 8.20 from (126). 

A rough lower bound on H(n) for any procedure can be easily obtained 
from information theory. The total information in n units is 

—n(p log 2 p + tf logo q), 

and this is to be equated with the product of the expected number of 
tests H{n) and the average information obtained per group-test. Since 
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the maximum information per group-test (in bits) is unity, we obtain 
for any procedure 

H(n) ^ -n(p log 2 p + q log 2 q). (127) 

For n = N = 100 and q = 0.99, this gives 8.09 as a lower bound. Since 
a result better than 8.09 is impossible, the smallness of the difference 
8.32 — 8.09 = 0.23 is an indication of how far Ri can possibly be from 
an optimal solution. However, it should not be inferred that the lower 
bound (127) can be reached for any value of q (except possibly for 
p = q = %) by any procedure. In fact, for q < \ and q decreasing to- 
wards zero, it has been shown that an optimal procedure must have 
H (to) = n, whereas the right member of (127) approaches zero. 

It has been pointed out to the authors by S. W. Roberts that a lower 
bound for G(m, n) for any procedure is easily shown to be 

+ (n - i) (p log 2 p + q log 2 q) 
= y^ W l log 2 q" + ( 1 - q m ) log 2 (1 - f)] {m) 
- in + "^ - J (p log 2 p + q log 2 q) . 

XIII. LACK OF OPTIMALITY OF PROCEDURE R x 

To illustrate the fact that Ri is not optimal in the general case when 
units are identifiable and "mixing" of units from the binomial and de- 
fective sets is allowed, we shall describe a method of obtaining an im- 
provement on R x . It is sufficient to consider the case N = 3, but the 
case N = 4 is more typical, and we shall use the latter. Let R * denote a 
procedure for N = 4, part of which is described by Fig. 2 and the remain- 
ing part of which is arbitrary. (We can therefore also regard R * as a 
set of procedures, with the common part shown in Fig. 2.) Let oi , a? , 
b { , bi denote individual units; it will be assumed that the a-units are 
distinguishable from the 6-units. The part of Fig. 2 enclosed by dashed 
lines is different from R\ , since it includes mixing; the rest of the pro- 
cedure agrees with R\ for q close to unity. For q close to unity and m = 2, 
after the first two group-tests result in failure, we should act as if there 
was exactly one defective present until it is proved otherwise. Then, for 
q close to unity, the above procedure R * terminates in one or two ad- 
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PROCEDURE R* (STARTING WITH H (4) -SITUATION) 



TEST 
NUMBER 



(a,,a 2f b,,b 2 ) 



END 




(ANY FINITE 
CONTINUATION) 



Fig. 2 — Procedure Ro*. 



ditional tests with probability close to \ for each. More precisely, we 
obtain for the conditional expected number of additional tests required 
under Ro*, given that the first two tests result in a failure, 



G(2,4|# *) = 



VQ 



1 - q< 



+ 



2pq 3 



+ 



P 2 f(q) 3g 3 + pf(q) 



1 - 5 2 



1 +<Z 



(129) 



where f(q) is a polynomial in q. In comparison, we have under Ri for 
0.707 < q < 1.000 



Gl(2 ,4) = l+^ + ^< 2 > 



1 - q 2 



1 - f 



1 



(130) 



1 + 5 



(6 + 2q - 2f/ - 2q a ). 



For q approaching unity, the value in (129) approaches f, while that in 
(130) approaches 2. This proves that any finite continuation in Fig. 2 
will be better than Ri for q sufficiently close to unity. In a particular 
procedure to be discussed in a separate paper, the dividing point for the 
(?(2, 4) -situation between "no-mixing" and "mixing" is 

q = (1 + V33)/8 = 0.843 

(to three decimal places). The maximum improvement over Ri for 
n = 4 in the expected number of tests required for the ^-situation is a 
decrease of 0.04. The price to be paid for this improvement will be an 
increase in the complexity of the procedure. 
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APPENDIX A 

The Information Procedure 

Another procedure which was investigated is based on choosing that 
value of x which maximizes the "amount of information" that the next 
test will give. The amount of information in a test with two outcomes is 
p log2( 1/p) + q log 2 ( \/q) if p is the probability of either outcome. Hence, 
equating the information in an //-situation obtained by taking x and 
x + 1 in the next test gives 

f/log 2 f/+ (1 - <niog 2 (l - q x ) 

= q x+1 log 2 q x+l + (1 - q x+l ) log 2 (1 - r/ +1 ), ^ 

whose root will be used as a dividing point between x and x + 1. It is 
easy to verify that, for any integer x and any integer n ^ x + 1, the 
unique positive root of (131) is also the unique positive root of 

1 - q x - q x+1 = 0. (132) 

It is interesting to note that the solution implicit in (132) is easily seen 
to be equivalent to finding that positive integer x (for fixed, known q) 
for which q x is closer to \ than is either q x ~ l or q x+1 . In fact, by (132), 
the right endpoint of the interval for x is such that q x and q x+l are 
centered about § and the left endpoint of the interval for x is such that 
q x and q x ~ are centered about \. Similarly, for the G-situation with 
m > 1, we equate 

(f^(t£) + (i^H(^) <») 

with the same expression, except that x is replaced by x + 1, and find 
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that, for any n ^ m > 1, the dividing point between x and x 4- 1 is 
the unique root in the interior of the unit interval (if it exists) of 

1 - Q x - q* +1 + <T = 0. (134) 

If we remove the root q = 1 in (134), the dividing point is the unique 
positive root (if it exists) of 

ii I 2 , x-l x+l x+2 m-l _ fx /lOK^ 

i + q + q • • ■ + q — q — q — • • • — q - u. uoo; 

If the root does not exist for some m > 1, then x + 1 will never be used 
for that m. It should be noted that the left member of (134) is a strictly 
increasing function of x and, for x ^ (m - l)/2, m > 1 and any 
fixed q with ^ g < 1, we have 

l-cf - j-h + ff - ^ (1 - <2 ( '- n/2 )(l - g (m+J)/z ) > 0. (136) 

It follows that the highest value of x for which a nondegenerate root ex- 
ists is such that x + 1 < (m + l)/2 and hence, under this procedure, 
we never take a test group of size greater than m/2. It is interesting to 
note that the dividing points for any G-situation do not depend on n. 
These equations define a new procedure R 2 , which we shall also call 
the information procedure. For this, let F 2 (m) denote the expected num- 
ber of group-tests required to "break up" a defective set of size m, i.e., 
to reach an tf-situation. Let F 2 *(m) = (1 - q'^p^Fzim). Then we 
can write as in (20), for any n and for the appropriate interval where the 
next test group is of size x, 

F 2 *(m | x) = E r/- 1 + q x F 2 *(m - x) + F 2 *(x) (m > 1), (137) 

H 2 ( n \ x) = 1+ qH 2 {n - x) + pF 2 *(x) + pT.q^lUn - i). (138) 

»=i 

The boundary conditions state that F 2 *(l) = H 2 (0) = for all q. Those 
expressions for F 2 *(m) which are used to generate expressions for H 2 (n) 
for 2 S m ^ n ^ 12 are given in Table VIII; the resulting expressions 
for H 2 (n) (with .r-values) are given in Table VI. Table VII gives the 
dividing points for n = 1(1)100 and for m = 1(1)16, 20(5)100 for 
procedure R 2 . 

It should be noted that the F 2 *(m) as well as the H 2 (n) -function are 
not all continuous. At the point of discontinuity the x corresponding to 
the smaller expectation should be used. 

It is interesting to observe in the numerical comparisons of Table III A 
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that the procedure R 2 compares quite favorably with the procedure Ri . 
In addition, the fact that the dividing points are easier to compute 
makes it better for practical applications, since the dividing points for 
R\ are only known exactly up to n = 16. It is also interesting to note 
that the limiting expressions in Table III A as n — > <=° and in Table III B 
asm-) co are the same as (132). 

It is interesting to note that a succession of modifications R2 3 (j = 1 , 
2, . . . ) of the information procedure, R 2 , are possible such that R2 1 = 
R 2 and Ri 3) = Ri for j ^ M(m,n). Here M(m,n), as defined in (10), is 
the maximum number of group-tests required if we start with a G(m,n)- 
situation [where G(0,n) corresponds to H(n)]. Under the procedure R 2 ' }) 
we find and use that x which maximizes the ratio of the information ex- 
pected from at most j group-tests to the conditional expected number 
of tests required given that we will stop after at most j group-tests. In 
the special case when there is no possibility of stopping before j tests, 
we can disregard the denominator and simply maximize the information. 
For the case j = 1 this is clearly equivalent to R2 . For j ^ M(m,n) the 
information expected from at most j tests is the same regardless of what 
x is used next and of what sample path is taken, since all units are then 
analyzed. Hence, the numerator above can be disregarded and the prob- 
lem is to minimize the denominator or expected number of tests. Under 
the assumption of "no-mixing" of units from the binomial and defective 
sets, this gives the procedure R\ . 

For any //(n)-situation with n ^ 4 and j ^ 2, these procedures appear 
to eliminate the possibility of taking n — 1 units in the next test-group. 
For example, if n — 4,j — 2 and q > 0.618, then we will want to com- 
pare x = 2 and x = 3. For j = 1, the dividing point between x = 2 
and x = 3 is q = 0.755. Since neither x = 2 nor x = 3 can result in 
termination after one test, we can disregard the denominator and com- 
pare for x = 2 and x = 3 the information expected from two group-tests. 
After simplification, the difference between the results expected after 
x = 2 and x = 3 can be written as 

tf-q[(l + q) log, (1 + q) - q log 2 q] ^ 0, (139) 

which shows that x = 2 is preferable to x = 3 for all q > 0.618. The 
same result holds for all j ^ 2. Then we find that the dividing point 
between x = 2 and x — 4 for j = 2 is the nondegenerate root between 
zero and unity of 

(2 _ 2 q - q 4 - r/ - 2q 7 )q log 2 q - (2 + q*) 

(1 - q") logs (1 + q) - q\l - q*) logs p = 0, 
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which is 0.789 to three decimal places. For j ^ 8 = M (0,4) the corre- 
sponding dividing point for R 2 ll) = Ri is the root of 1 — q — q = 
0, or 0.786 to three decimal places. Curiously enough, the same result 
0.786 is also the dividing point between x = 2 and x = 4 fori = 1. 

For the special case x = 1, we state without proof that, for any 
H(n)-situation, the dividing point between x — 1 and x = 2 is again 
KV5 — 1) = 0.618 to three decimal places. 

Formulae for the expected number of tests under #2 3 f or 1 < j < 
M(m,n) have not been derived in this paper. 

APPENDIX B 

Definition of Procedure R s 

It may happen in some problems that recombination is undesirable or 
impossible, or it may be that we are interested in finding out just how 
much is saved by allowing recombinations. Both are good reasons for 
considering a procedure R3 that is similar to Ri except that recombina- 
tions are not allowed. This simply means that any two operationally 
formed sets cannot be combined to form a new set from which subse- 
quent test-groups are to be taken. In procedure Ri the possibility of 
mixing proper subsets of two different sets was never used, and the same 
will be true for R 3 . If both recombinations and mixing are not used, 
then, as the experiment continues, the operationally formed sets can 
only be broken down further into smaller and smaller sets, yielding a 
nested set of partitions; i.e., any two units separated at some stage re- 
main separated in subsequent stages. It follows that any defective set 
present is not affected by the number or size or nature of other sets 
present. Hence, we define G 3 (m), with a single argument, as the condi- 
tional expected number of group-tests required to remove all the de- 
fectives from a set of size m which is known to have at least one defec- 
tive. 

The recursion formulae for R 3 are, for n 2: 1 and m > 1, 

H 3 (n) = 1 + min {q x H 3 (n — x) 



1 < J- < n 



(141) 
+ (1 -q*)[G*(x) +H 3 (n-x)}}, 

G,(m) = 1 + min |Y f - Q G 3 (m - x) 

iM-^ IVl " 5 ) ^ (M2) 
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with boundary conditions H 3 (0) = (? 3 ( 1 ) = for all q. If we let G 3 *{m) 
denote (1 — q m )p~ G$(m) and simplify, we obtain 

H 3 (n) = 1 + min \H 3 (n - x) + pG 3 *(x)}, (143) 

lgxgn 

m r 

G 3 *(m) = Z <T X + min {(?,*(*) + q x G 3 *(m - x) 

+ H 3 (n - x) (t, g" 



l<x<m— 1 



(144) 



with boundary conditions H 3 (0) = G 3 *(l) = for all q. 

Numerical comparisons of the results for R t and R 3 are given in Ta- 
bles II A and II B. 



APPENDIX C 



Definition of Two Halving Procedures 

Two "halving" procedures R A and R 5 are defined below, the principal 
purpose being to compare the results on the expected number of tests 
required with comparable results for Ri and R 3 . The procedure Ri al- 
lows recombinations exactly as in Ri , while R b is the same as R4 except 
that recombinations are not allowed. Both R4 and R$ are of particular 
interest, since they can be carried out without knowing the true value 
of q. 

The procedure R A is carried out like Ri except that, if the defective 
set is of size m > 1, the next test group is a subset of size m' = [m/2] 
(i.e., the largest integer contained in m/2) randomly selected from the 
defective set and, if m = or 1, the entire binomial set is used in the next 
test-group. In particular, we start with all iV units in the first test-group. 
The recursion formulae for J? 4 are 

Ht(n) =g n -r-(l-g n )[l+G 4 (n,n)] = l-f-(l - <f)G 4 (n,n), (145) 

(? 4 (m,n) = 1 + ( \ - q m ) Gt(m, n) 

m , (1 _ q —'\ , , (146) 

+ Q m ( — : — ) Gi(m — m ,n — m) (m > 1 ), 

with the same boundary conditions as in Ri . If we let Fi(m) denote the 
expected number of tests required to break up a defective set of size m, 
it can be shown as in the case of R\ that 
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(? 4 (m, n) = F<(m) + (r^-^) E « W ff«(n - 0- (147) 
\i — g / i=i 

If we let 

G 4 *(m,n) -(^=£V(m J n) and /?**(«) =(V=^") F * (m) ' (148) 

the recursion formulae for i? 4 reduce to 

F*(m) = £ g'- 1 + F 4 *(m') + g m F 4 *(ro - m') (m > 1), (149) 

ff 4 (n) = 1 + pG 4 *(m, n) 

(150) 
= 1 + p{F 4 *(n) +D q^Htin - *)}, 

with boundary conditions F 4 *(l) = # 4 (0) = for all q. For n g 5, 
the results are the same as those for #i , if we take q close to unity 
in the formulae for fi L . The results for H 4 (n) for n = 6(1)12 are, for 
all?: 

tf 4 (6) = 14 - 9? - 2 9 2 - g 3 - g 5 , 

H 4 (7) = 17 - llg - 3g 2 - q + g 4 - 2? 5 , 

tf 4 (8) = 21- 14g-4g 2 -g 3 4-g 4 -2g 5 , (151) 

ff 4 (9) = 25 - 18g - 4g 2 - q + q - 2g 5 + q - q\ 

Hi(10) = 29 - 22g - 4? 2 + q - 3g 6 + q - q, 

ff 4 (ll) = 33 - 26tf - 4 ? 2 + <1 ~ 4( f + V + (2 7 - S 8 - ? 9 > 

ff 4 (12) = 37 - 29g - 4g 2 - 4r/ + 2 5 6 + q - q % - q. 

The recursion formulae for the halving procedure R b , which continues 
to separate sets into smaller and smaller subdivision, are 

ff.(n) = 1 + (1 - q n )G 6 (n), (152) 



(153) 



0.(m) = i + (rrr^r) [G&(m/) + Hi{m ~ m ' )] 

(. m—m'\ 

1 ~_ g m j G 6 (m - w') (m>l), 

where ra' is defined as above. Here it was not necessary to use a double 
argument with G because there is no recombination allowed. The bound- 
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ary condition is G&(1) = for all q. If we define G*{m) as in (153) 
and use (152), the recursion formulae reduce to 

H 6 (n) = 1 + pG b *(n), (154) 



G B *(m) =E? M +E9 H 



+ G B *(w / ) + G b *(m - m') (m > 1 



(i5o; 



with the boundary condition G*{\) = for all q. For n ^ 3, the results 
are the same as for Ri if we take q close to unity in the formulae for Ri . 
The results for H 6 (n) f or n = 4(1)12 are: 



ff.(4) = 7 - 2g - 3g 2 - q\ 








H 6 (5) = 9 - 3g - 3g 2 - g 3 


- q\ 






H B (6) = 11 - Aq - 2g 2 - 3? 3 - 




6 

- Q , 




H b (7) = 13 - 4g - 4f/ - 2 S 3 - r/ 




- q\ 




H B (8) = 15 - 4g - 6g 2 - 3g 4 






8 

- q, 


// 5 (9) = 17 - 5g - 6g 2 - g 3 - 2q - 


- q 




9 

- q, 


ff 6 (10) = 19 - 6g - 6g 2 - 2g 3 


-3g 5 




10 

- q : 


/7 6 (H) = 21 - 7g - 5g 2 - 4g ; ' 


2g 6 - 


n 

- q 


11 

- q : 


/7 5 (12) = 23 - 8g - 4g 2 - 6g 3 




- 3g 


12 

- q ■ 



:i56) 



APPENDIX D 

Known Procedures 

An attempt has been made to put the Dorfman procedure and the 
Sterrctt procedure 3 in the best form that is comparable with the other 
procedures treated here. For each of N = 4, 8 and 12, we have found 
the division into equal (or approximately equal) subsets such that 
the Dorfman plan of testing defective sets one at a time gives the small- 
est possible expected number of tests required. It should not be inferred 
that these results would be the same if a straightforward application of 
the tables published by Dorfman and Sterrctt, respectively, were made, 
since their tables are only concerned with very large iV. In the Dorfman 
plan we use a common test-group size for binomial sets and, for defective 
sets, the units are all tested one at a time. In the Sterrett plan, there is a 
common test group size for binomial sets at the outset and, for defective 
sets, the units are tested one at a time only until a defective unit is found. 
Then the remaining units, from that defective set only, are pooled and 
tested. This is continued until that particular defective set is completely 
analyzed before we start with other sets. We have also assumed that 
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logical inference would be used whenever possible in the Dorfman and 
Sterrett procedures. 

For the Dorfman procedure R 7 , if the common group size is c, then, 
for any binomial set of size n (where n ^ c), we obtain 

H,M -f + l i-,)(n + l- j^l) (157) 

= n + 1 - g n_1 - (n - l)q n . 

For example, for N = 12 and g = 0.90, we find that d = c 2 = c 3 = 4 
gives the best results and, using (157), with c = 4, we obtain 

# 7 (12) = 3i/ 7 (4) = 3(5 - q 3 - 3g 4 ) = 6.908. (158) 

For the Sterrett procedure R 6 , if the group size at the outset is c, 
then, for each binomial set of size n (where n ^ c), we obtain 

H 6 (n) = q" + p[2 + H 6 (n - 1)] 4- gp[3 4- tf 6 (n - 2)] 

+ ••• + q n ~ 2 p[n + ff,(l)] 4- <T> 

n _ "-1 ._ (159) 

= g" + wpfy" ' + pE *V 2 + P Z ?' '^(n - *), 



i=2 «'=1 



with boundary condition H 6 (l) = 1 for all q. It can be verified (we omit 
the details) that the solution of this system is given by 

H 6 (n) = (2n - 1) - (n - l)g - E 3*'- (160) 

1=2 

For example, for N = 12 and # = 0.90, we find that c x = c 2 = c 3 = 4 
gives the best results and, using (160), with c = 4, we obtain 

#e(12) = 3ff«(4) = 6.315. (161) 

APPENDIX E 

Cost Considerations 

In this Appendix we introduce another procedure R s , which brings 
into play the cost of throwing away a good unit and balances it against 
the cost of conducting another group-test. It is interesting to note that 
R s was the solution given when the problem was first brought to the 
authors' attention in a practical application. 

For procedure R s we divide all the N units into approximately equal 
subsets of size x', where x' is the nearest positive integer to the solution 
in x of 
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(1 " VY = *, (102) 

where p is the known a priori probability of a unit being defective. It 
is assumed here that p « £. Each subgroup is tested either simultane- 
ously or in sequence, and good subgroups are removed. Assuming x' > 1, 
subsets shown to contain at least one defective are pooled. Since the 
same number of defectives has now been put into a pooled set of approxi- 
mate size N/2, it follows that the probability of drawing a defective from 
the pooled set is approximately double the original a priori probability 
p. Then the pooled subset is again divided into approximately equal sub- 
sets of size x", where x" is the nearest positive integer to the solution of 

(1 - 2p) x = i (163) 

The process is repeated (say) a total of t times. If tp gets larger than \, 
x is taken to be unity and the units would then be tested one at a time. 
However, it may be more economical to stop the procedure before tp 
reaches \ and scrap the pooled defective subgroup. The amount of sav- 
ing may be substantial if the cost of manufacturing a unit, say Co , is 
small compared to the cost of each group-test, say C\ . 

Suppose, for example, that we start with N = 8 units and q is given 
to be 0.90. The approximate solution of (162) is x' = 7 but, since this 
leaves a subset of size 1 , we make our first test on all 8 units. If the test 
is a success we are through; otherwise we look for a solution of (163) 
and find that x" = 3. The 8 units are divided into subsets of size 3, 3 
and 2, and each is tested. Good subsets are removed. We purposely 
avoid the next stage which requires testing one at a time. Hence, any 
of these three sets of size 3, 3, and 2 that proves to contain at least one 
defective is scrapped. 

Let T denote the number of tests required in the above example and 
let D denote the total number of scrapped units (i.e., good units and de- 
fectives that are discarded). Then, for this procedure R 8 with N = 8 
and q = 0.90, we obtain 

E\T | Rs\ = 4 - 3r/ = 2.709, (164) 

E\D | R s \ = 8 - 2q- - Gq 3 = 2.006. (165) 

If we define the expected loss E\L \ R,} for any procedure /?,■ by 

E\L | /?,) = c E\D | R t ) + ciE{T \ /?,}, (166) 

we find for procedures R } and R s , respectively, 

E\L | /?,( = 0.800 co + 3.904 c, , (167) 

E\L | 7? 8 1 = 2.006 co + 2.709 c, . (168) 
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A comparison of these two expressions shows that R s will be more eco- 
nomical in this case if the ratio 

Cj 2.006 - 0.800 = = j approximately. (169) 

Co _ 3.904 - 2.709 

Hence, R» is more economical in this case if the total cost of a single test 
is greater than the total cost of manufacturing a single unit. 

Similarly, it can be shown that it would be more economical to stop 
after the first test (and scrap all 8 units if there is at least one defective 
present) when the ratio of the two costs in (169) is greater than 1.492 
(or approximately 1.5). 
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Fig. 3 — Diagram showing number of observations to be taken in any H -situa- 
tion for n = 1 through 16 — for procedure lt\ . 
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Fig. 4 — Diagram showing number of observations to be taken in any G-situa- 
tion for to = 2 through 16 and any n ^ m — for procedure R i . 
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Table II A — Comparison of the Expected Number of Group-Tests 

Required for Different Procedures for N = 4, 8, 12 

and Selected Values of q 





Q 




0.00 


0.50 


0.75 


0.90 


0.95 


0.98 


0.99 


1.00 


(Proposed Procedure) 


4.000 
8.000 

12.000 


4.000 
8.000 
12.000 


3.332 
6.619 
9.905 


2.051 
3.904 
5.790 


1.538 
2.499 
3.594 


1.218 

1.612 
2.073 


1.110 

1.308 
1.543 


1.000 
1.000 
1.000 


R-, 
(Information Procedure) 


4.000 
8.000 
12.000 


4.000 
8.000 
12.000 


3.332 
6.663 
9.956 


2.051 
4.141 
5.839 


1.538 
2.500 
3.599 


1.218 
1.612 
2.078 


1.110 
1.308 
1.544 


1.000 
1.000 
1.000 


(Proposed without Re- 
combinations) 


4.000 

8.000 

12.000 


4.000 
8.000 
12.000 


3.375 

6.719 

10.062 


2.105 
4.052 
6.103 


1.576 
2.631 
3.851 


1.236 
1.680 

2.2011 


1.119 
1.345 
1.617 


1.000 
1.000 
1.000 


Ra 

(Halving with Recom- 
binations) 


4.000 
21.000 
37.000 


4.000 
12.S75 
21.408 


3.332 

7.670 

12.365 


2.051 
3.906 
6.021 


1.53S 
2.500 
3.620 


1.218 
1.612 
2.078 


1.110 

1.308 
1.544 


1.000 
1.000 
1.000 


R, 
(Halving without Recom- 
binations) 


7.000 
15.000 
23.000 


5.1S8 
11.309 
17.203 


3.496 

7.576 

11.653 


2.114 
4.141 
6.309 


1.578 
2.678 
3.900 


1.236 
1 . 700 
2.229 


1.119 
1.355 
1.627 


1.000 
1.000 
1.000 


R 6 
(Sterrett) 


4.000 

8.000 

12.000 


4.000 
8.000 
12.000 


3.375 

6.750 

10.125 


2.105 
4.210 
6.315 


1.576 
3.151 
4.333 


1.236 
1.807 
2.973 


1.119 
1.412 

1.852 


1.000 
1.000 
1.000 


Ri 

(Dorfman) 


1.(1(11) 

8.000 

12.000 


4.000 
8.000 
12.000 


3.375 

6.750 

10.125 


2.303 
4.605 
6. 90S 


1.699 
3.398 
5.097 


1.292 
2.176 
3.334 


1.148 
1.609 
2.354 


1.000 
1.000 
1.000 



Table II B — Comparison of Hi(n) and Hz(n) for Three Values 
of q and Larger n- Values 





9 = 


0.90 


i = 


0.95 


Q = 


0.99 




Hi(n) 


ffl(») 


Bi{n) 


B,(n) 


Bi(n) 


n,(„) 


10 


4.872 


5.101 


3.039 


3.242 


1.425 


1.481 


20 


9.572 


10.155 


5.940 


6.456 


2.051 


2.221 


30 


14.301 


15.209 


S.791 


9.626 


2.738 


3.057 


40 


19.024 


20.260 


11.671 


12.798 


3.478 


3.943 


50 


23.750 


25.361 


14.555 


16.009 


4.243 


4.853 


60 


28.475 


30.415 


17.438 


19.246 


5.026 


5.792 


70 


33.200 


35.469 


20.316 


22.415 


5.830 


6.754 


80 


37.925 


40.520 


23.197 


25.591 


6.647 


7.717 


90 


42.650 


45.621 


26.078 


28.781 


7.477 


8.683 


100 


47.375 


50.675 


28.959 


32.019 


8.230 


9.687 
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Table IV A — Formulae for Fi*(m) for Procedure Ri and Values 
of the Next Test-Group Size x for m = 2(1)16 





g-interval 


.V 


l 
1 


1 


•r- 


V 3 


5 4 


<7 5 


<i" 


Q 1 


gB 


<7» 


5 io 


<7 11 


«w 


git 


g» 


ff lS 

t 


Fi*(2) 


0.000 


to 1.000 


1 


fi*(3) 


0.000 


to 1.000 


1 


1 


2 


2 




















— 


F,*(4) 


0.000 


to 0.755 


1 


1 


2 


3 


3 


















0.755 


to 1.000 


2 


2 


2 


2 


2 






















— 


""I 


*V(5) 


0.000 


to 0.682 


.1 


1 


2 


3 


4 


4 
















0.682 


to 1.000 


2 


2 


2 


2 


3 


3 




















— 


"" t 




0.000 


to 0.652 


1 


1 


2 


3 


4 


5 


5 














fi*(8) 


0.652 


to 0.755 


2 


2 


2 


2 


3 


4 


4 




















t 




0.755 


to 1.000 


2 


2 


2 


3 


3 


3 


3 


















— 


— — 




0.000 


to 0.637 


] 


1 


2 


3 


4 


5 


6 


6 














W7) 


0.637 


to 0.682 


2 


2 


2 


2 


3 


4 


b 


b 




















0.682 


to 0.857 


2 


2 


2 


3 


3 


3 


4 


4 






















0.857 


to 1.000 


3 


2 


3 


3 


3 


3 


3 


3 

















— 


"I 




0.000 


to 0.629 


1 


1 


2 


3 


4 


5 


6 


7 


7 














0.629 


to 0.652 


2 


2 


2 


2 


3 


4 


5 


6 


(i 
















t 


^i*(8) 


0.652 


to 0.755 


2 


2 


2 


3 


3 


3 


4 


5 


5 
















t 


0.755 


to 0.819 


2 


2 


2 


3 


3 


4 


4 


4 


4 
















t 




0.819 


to 0.899- 


3 


2 


3 


3 


3 


3 


3 


4 


4 




















0.899" 


to 1.000 


4 


3 


3 


3 


3 


3 


3 


3 


3 














— 


t 




0.000 


to 0.624 


1 


1 


2 


3 


4 


5 


6 


7 


8 


8 












0.624 


to 0.637 


2 


2 


2 


2 


3 


4 


5 


6 


7 


7 














t 


Fi*(9) 


0.637 


to 0.682 


2 


2 


2 


3 


3 


3 


4 


5 


6 


6 














t 


0.682 


to 0.S03 


2 


2 


2 


3 


3 


4 


4 


4 


5 


5 














t 




0.803 


to 0.881 


3 


2 


3 


3 


3 


3 


4 


4 


4 


4 














t 




0.881 


to 1.000 


4 


3 


3 


3 


3 


3 


3 


3 


4 


4 


















0.000 


to 0.622 


1 


1 


2 


3 


4 


5 


6 


7 


8 


9 


9 












■• 




0.622 


to 0.629 


2 


2 


2 


2 


3 


4 


5 


6 


7 


8 


8 












T 




0.629 


to 0.652 


2 


2 


2 


3 


3 


3 


4 


5 


6 


'/ 


7 














fi*(10) 


0.652 


to 0.755 


2 


2 


2 


3 


3 


4 


4 


4 


5 


6 


6 












t 




0.755 


to 0.782 


2 


2 


2 


3 


3 


4 


4 


b 


5 


b 


b 
















0.782 


to 0.857 


3 


2 


3 


3 


3 


3 


4 


4 


4 


b 


5 












t 




0.857 


to 1.000 


4 


3 
1 


3 

2 


3 
3 


3 

4 


3 
5 


3 

6 


4 

7 


4 

s 


4 
9 


4 
10 


10 


— 






— 


"" t 




0.000 


to 0.620 


1 




0.620 


to 0.624 


2 


2 


2 


2 


3 


4 


5 


6 


V 


8 


y 


9 










t 




0.624 


to 0.637 


2 


2 


2 


3 


3 


3 


4 


b 


6 


7 


8 


8 










t 


fi*(ll) 


0.637 


to 0.6S2 


2 


2 


2 


3 


3 


4 


4 


4 


5 


6 


7 


7 










t 


. 682 


to 0.779 


2 


2 


2 


3 


3 


4 


4 


5 


5 


5 


6 


6 










\ 




0.779 


to 0.819 


3 


2 


3 


3 


3 


3 


4 


4 


5 


5 


5 


b 












0.819 


to 0.857 


3 


2 


3 


3 


3 


4 


4 


4 


4 


4 


5 


5 










t 




0.857 


to 1.000 


4 


3 


3 


3 


3 


3 


4 


4 


4 


4 


4 


4 








— 


t 




0.000 


to 0.619 


1 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


11 








0.619 


to 0.622 


2 


2 


2 


2 


3 


4 


5 


6 


7 


8 


9 


10 


10 








t 




0.622 


to 0.629 


2 


2 


2 


3 


3 


3 


4 


5 


6 


7 


8 


y 


9 








t 




0.629 


to 0.652 


2 


2 


2 


3 


3 


4 


4 


4 


5 


6 


7 


8 


8 








t 


*\*(12) 


0.652 


to 0.755 


2 


2 


2 


3 


3 


4 


4 


5 





b 


6 


7 


7 








t 


0.755 


to 0.768 


2 


2 


2 


3 


3 


4 


4 


5 


5 





6 


6 


6 








' 




0.768 


to 0.803 


3 


2 


3 


3 


3 


3 


4 


4 


5 


5 


5 


6 


6 








t 




0.803 


to 0.846+ 


3 


2 


3 


3 


3 


4 


4 


4 


4 


5 


b 


b 


5 








t 




0.846+ to 0.899- 


4 


3 


3 


3 


3 


3 


4 


4 


4 


4 


4 


b 


5 












0.899- 


to 1.000 


4 


3 


3 


3 


3 


4 


4 


4 


4 


4 


4 


4 


4 











t These equations are not needed to compute H x {n) -formulae if the experiment 
starts with a "pure binomial" set of any size N (i.e., if m = at the outset). 





Table 


T\ 


r A- 


- Continued 


















(/-interval 


X 

1 


l 
1 


2 


3 


■/ 3 
4 


5 


9* 

() 


,6 

7 


s 


9 


10 


•I'- 
ll 
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The exponential symbols -f- , — indicate only the relative magnitude of two 
different roots that are equal to three decimal places (i.e., ar < a+). 

1231 



1232 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1959 



Table IV B — Formulae for Hi(n) and Values of the Next Sample 
Size x for Any //-situation Expressed as Functions of 
q for n = 2(1)12 When the Procedure Ri Is Used. 
The integer shown below q" opposite H\(n) is the coefficient of q" in the poly- 
nomial formula for H\(n). 
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to 0.819 


3 


26 


-17 


-7 


2 


2 


-5 


2 


2 


-2 





1 








0.819 


to 0.824 


4 


30 


-24 


-4 





4 


-3 


-1 


1 


1 





-1 








0.824 


to 0.840 


4 


31 


-26 


-3 





4 


-3 


-1 


1 


1 


-1 










0.840 


to 0.844 


4 


31 


-26 


-3 


1 


3 


-4 





1 


1 


-1 


-1 


1 




#i(12) 


0.844 


to 0.856 


4 


32 


-28 


-1 





2 


-3 


-1 


2 


1 


-1 


-1 


1 




0.856 


to 0.857 


4 


32 


-28 


-1 


1 





-2 


-1 


3 





-2 








1 




0.857 


to 0.869 


4 


32 


-28 


-1 


2 


-2 


-1 


-1 


2 


1 


-1 


-1 





1 




0.869 


to 0.881 


6 


32 


-27 


-1 


2 


-3 


-2 





1 


2 


-1 


-1 








0.881 


to 0.885 


6 


32 


-27 








-2 


-2 








3 





-2 








0.885 


to 0.896 


6 


32 


-26 


-1 





-3 


— 1 


1 


-1 


1 


2 


-2 








0.896 


to 0.899- 


6 


33 


-27 


-1 


-1 


-2 


— 1 


1 
















0.899- 


to 0.904 


6 


34 


-29 


1 


-3 


-1 


— 1 





1 








1 


-1 






0.904 


to 0.912- 


6 


35 


-30 





-2 


-1 


— 1 


1 











-1 


1 






0.912- 


to 0.912+ 


7 


35 


-29 


-1 


-2 


-1 


— 1 


1 





-1 








1 






0.912 + 


to 0.919 


7 


36 


-30 


-2 


-1 


-1 











-1 


-1 


1 





1 




0.919 


to 1.000 


12 


37 


-31 


-2 


-1 














-1 


-1 









The exponential symbols +, — indicate only the relative magnitude of two 
different roots that are equal to three decimal places (i.e., a~ < a + ). 



Table IV C — Formulae for Gi*(m, n) and the Values of the 

Next Sample Size x for Certain G-Situations Expressed as 

Functions of q for n = 2(1)12 and q ^ 0.850 When 

Procedure Ri is Used 

The integer shown below q v opposite Oi'(m,n) is the coefficient of q y in the polynomial formula for 
<7i*(m,m). 



n 


m 


g-interval 


.T 

1 


(7i*(m,ff) 


1 

2 


'/ 


<l" 


«' 


<7 4 


? 5 


9° 


9 7 


9* 


<7° 


qia 


8" 


2 


2 


0.850 to 1.000 


1 






















3 


2 
3 


0.850 to 1.000 
0.850 to 1.000 


1 

1 


4 
4 


1 
2 


-1 
1 




















4 


2 
4 


0.850 to 1.000 
0.850 to 1.000 


1 

2 


6 

7 


2 
3 


-2 
1 


-2 


















5 


2 

: 


0.850 to 1.000 
0.850 to 1.000 
0.850 to 1.000 


1 
1 

2 


9 

9 

10 


2 

3 
3 


-4 
1 
1 


_2 

-3 

1 


-1 

-2 

1 
















6 


2 
3 
6 


0.850 to 1.000 
0.850 to 1.000 
0.850 to 1.000 


1 
1 
2 


12 
12 
13 


2 
3 
3 


-6 

1 
2 


-2 
-4 

2 


-1 

-2 
1 


-1 

-2 














7 


2 
3 

4 

7 


0.850 to 1.000 
0.850 to 1.000 
0.850 to 1.000 

0.850 to 0.857 
0.857 to 1.000 


1 

1 
2 

2 
3 


15 
15 
10 

10 
10 


2 
3 
3 

3 
4 


-8 
2 
2 

3 
3 


-2 

-6 

1 

2 

2 


-1 

-3 
-5 

1 
1 


-1 

-2 
-3 

1 


-1 

-1 
-2 

1 












8 


2 
3 
4 

8 


0.850 to 0.869 
0.869 to 1.000 

0.850 to 0.869 
0.869 to 1.000 

0.850 to 0.869 
0.869 to 1.000 

0.850 to 0.869 
0.869 to 0.899" 
0.899" to 1.000 


1 
1 

1 
1 

2 

2 

3 
3 
1 


18 

is 

18 
18 

111 
19 

19 
19 
20 


2 
3 

3 
4 

3 

1 

•1 
5 
5 


-11 
-11 

2 
2 

2 
2 

3 
3 
3 


-2 
_2 

-9 
-9 

1 
1 

2 
2 
2 


-1 
-1 

-3 
-3 

-7 
-7 

1 
1 
1 


-1 
-2 

-1 
-2 

-3 
-4 

1 



-1 
-1 

-1 
-1 

-2 
-2 

1 
1 


1 

-1 

-1 

2 

1 













2 
3 
•1 


0.850 to 0.869 
0.869 to 0.885 
0.885 to 0.899" 
0.899- to 1.000 

0.850 to 0.869 
0.869 to 0.885 
0.885 to 0.899" 
0.899- to 1.000 

0.850 to 0.869 
0.869 to 0.885 
0.885 to 0.899- 
D.899- to 1.000 


2 
2 
2 
2 


21 
21 
21 
22 

21 
21 

21 
22 

22 
22 
22 
23 


2 
3 
1 
3 

3 
4 

5 
4 

3 
4 

■ r > 
1 


-14 

-14 
-14 
-14 

2 

2 
2 
2 

2 
2 
2 

2 


-2 
-2 
-2 
-2 

-12 
-12 
-12 
-12 

1 

1 
1 

1 


-1 
-1 

-2 
-2 

-2 
-2 
-3 
-3 

-9 

-9 

-10 

-10 


-1 
-2 
-2 
-2 

-1 
-2 

-2 
-2 

-3 
-4 
-4 
-4 


-1 

-1 



-1 

-2 
-2 
-1 
-2 

-2 
-2 
-1 
-2 


1 





-1 
-1 

-1 

-1 
-1 
-1 


1 

1 

-1 

1 

1 

-1 

-2 
-1 
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Table I\ 


? c- 


- Continued 
















q-inUrval 




Gi'Un.n) 


« 


HI 


.V 


































1 


1 


<7 2 


9 3 


g* 


<7 6 


«« 


?' 


q* 


? 9 


g 10 


ff" 


11 




J. 850 toO. 809 


2 


22 


3 


2 


2 


2 


-7 


-4 


-1 












■ 


5.869 to 0.885 


■> 


22 


1 


2 


2 


2 


-8 


-4 


-2 












0.885 to 0.899" 


■> 


22 


5 


2 


2 


1 


-8 


-3 


-2 


-2 












J.899- to 1.000 


i 


23 


4 


2 


2 


1 


-8 


-4 


-2 


-1 












3.850 to 0.869 


3 


22 


4 


3 


2 


2 


2 


1 


2 


2 












0.869 to 0.881 




22 


5 


3 


2 


2 


1 


1 


1 


2 










9 


0.881 to 0.885 




23 


5 


3 


2 


2 








1 


2 












0.885 to 0.899" 




23 





3 


2 


1 





1 


1 














0.899- to 1.000 




24 


5 


3 


2 


1 


() 


» 


1 


1 








10 


0.850 to 0.857 




24 


2 


-17 


—2 





-l 


-1 


1 


1 


-i 










0.857 toO. 809 




24 


2 


-17 


-1 


-1 


-l 


-1 





1 










2 


0.869 to 0.885 




24 


3 


-17 


-2 


-1 


-2 


-1 





1 


i 








0.885 to 0.896 




24 


1 


-17 


-2 


-2 


-2 








-1 


i 










0.890 to 0.899- 




25 


1 


-17 


-3 


-2 


-2 





1 


-1 


-l 










0.899- to 1.000 




26 


3 


-17 


-3 


-2 


_2 


-1 


1 


(1 


-i 










0.850 to 0.857 




24 


3 


2 


-15 


-1 


-2 


-2 


1 


1 












0.857 to 0.809 




24 


3 


2 


-14 


-2 


-2 


-2 





1 


i 








3 


0.809 to 0.885 




24 


4 


2 


-14 


-2 


-3 


-2 


-1 


1 


i 








0.885 to 0.896 




24 


5 


2 


-14 


-3 


-3 


-1 


-1 


-1 


i 










0.896 to 0.899- 




25 


5 


2 


-15 


-3 


-3 


-1 





-1 


-l 










0.899" to 1.000 




20 


4 


2 


-15 


-3 


-3 


-2 








-l 










0.850 to 0.857 


2 


25 


3 


2 


1 


-11 


-3 


-2 
















0.857 to 0.869 


2 


25 


3 


2 


2 


-12 


-3 


-2 


-1 





i 








4 


0.869 to 0.885 


2 


25 


4 


2 


2 


-12 


-4 


-2 


-2 





i 








0.885 to 0.896 


2 


25 


5 


2 


2 


-13 


-4 


-1 


-2 


-2 


i 










0.896 to 0.899- 


2 


26 


5 


2 


1 


-13 


-4 


-1 


-1 


-2 


-l 








0.899" to 1.000 


2 


27 


4 


2 


1 


-13 


-4 


-2 


-1 


-1 


-i 










0.850 to 0.857 


2 


25 


3 


2 


2 


;i 


-10 


-4 








-l 










0.857 to 0.869 


2 


25 


3 


2 


3 


2 


-10 


-4 


-1 












5 


0.869 to 0.885 


2 


25 


4 


2 


3 


2 


-11 


-4 


-2 












0.885 to 0.896 


2 


25 


5 


2 


3 


1 


-11 


-3 


-2 


-2 












0.896 to 0.899" 


2 


26 


5 


2 


2 


1 


-11 


-3 


-1 


-2 


-2 










0.899" to 1.000 


2 


27 


4 


2 


2 


1 


-11 


-4 


-1 


-1 


-2 










0.850 to 0.857 


2 


26 


3 


3 


2 


3 


1 


-8 


-2 


-1 


-1 










0.857 to 0.809 


2 


211 


3 


3 


3 


2 


1 


-8 


-3 


-1 












O.809 to 0.885 


2 


25 


4 


3 


3 


2 





-8 


-4 


-1 










6 


0.885 to 0. 890 


2 


25 


5 


3 


3 


1 





-7 


-4 


-3 












0.890 to 0.899" 


2 


26 


5 


3 


2 


1 





-7 


-3 


-3 


-2 










0.899- to 1.000 


2 


27 


4 


3 


2 


1 





-8 


-3 


-2 


-2 










0.850 to 0.857 


3 


25 


4 


3 


2 


3 


2 


1 


3 


3 


2 










0.857 to 0.869 




26 


4 


3 


3 


2 


1 


1 


2 


2 


2 








10 


0.869 to 0.885 




26 


5 


3 


3 


2 





1 


1 


2 


2 








0.885 to 0.890 




20 


6 


3 


3 


1 





2 


1 





2 










0.896 to 0.899- 




27 





3 


2 


1 





2 


2 














0.899" to 1.000 




28 


5 


3 


2 


1 





1 


2 


1 








11 




0.850 to 0.857 




27 


2 


-20 


-1 





-1 


-2 


1 


2 


-1 


-1 








0.857 to 0.869 




27 


2 


-20 





-1 


-1 


-2 





2 





— 1 








0.869 to 0.885 




27 


3 


-20 


-1 


-1 


-2 


-2 





2 




- 1 






2 


0.885 to 0.890 




27 


4 


-20 


-2 


-2 


-2 








-1 












0.890 to 0.899" 




28 


4 


-20 


-3 


-2 


-2 





1 


-1 


-1 










0.899" to 0.904 




29 


3 


—20 


-3 


-2 


-2 


-1 


1 





-1 










0.904 to 1.000 




30 


3 


-21 


-3 


-2 


-2 





1 





1 _1 
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Table IV C - 


— Continued 


























Gi* 


m,n) 










ft 


m 


(/-interval 


X 






















i 

27 


1 


<f- 


q> 


g* 


q* 


<7 6 


9 7 


Q* 


9 9 


„io 


</" 


11 




0.850 to 0.857 


3 


2 


-17 


-l 


—2 


-2 


1 


1 





-1 








0.857 to 0.869 




27 


3 


2 


-16 


_2 


-2 


-2 





1 


1 


-1 








0.869 to 0.885 




27 


4 


2 


-16 


-3 


-3 


_2 


-1 


2 


1 










0.885 to 0.896 




27 


5 


2 


-16 


-4 


-3 


-1 


-1 





1 








3 


0.896 to 0.899" 




28 


5 


2 


-17 


-4 


-3 


-1 








-1 










0.899' to 0.904 




29 


4 


3 


-18 


-4 


-3 


—2 








-1 


1 








0.904 to 1.000 




30 


4 


2 


-18 


-4 


-3 


-1 








-1 


-1 








0.850 to 0.857 


2 


28 


3 


2 


2 


-14 


-3 


-3 

















0.857 to 0.869 


2 


28 


3 


2 


3 


-15 


-3 


-3 


-1 




1 










0.869 to 0.885 


2 


28 


4 


2 


3 


-15 


-4 


-3 


-2 




1 








■1 


0.885 to 0.896 


2 


28 


5 


2 


3 


-16 


-4 


-2 


-2 


—1 


1 










0.896 to 0.899" 


2 


29 


5 


2 


2 


-16 


-4 


_2 


-1 


_1 


-1 










0.899- to 0.904 


2 


30 


4 


3 


1 


-16 


-4 


-3 


-1 


_1 


-1 


1 








0.904 to 1.000 


2 


31 


I 


2 


1 


-16 


-4 


-2 


-1 


-1 


-1 


-1 








0.850 to 0.857 


2 


28 


3 


2 


3 


3 


-13 


-4 








-1 










0.857 to 0.869 


2 


28 


3 


2 


4 


2 


-13 


-4 


-1 














0.869 to 0.885 


2 


28 


4 


2 


4 


2 


-14 


-4 


-2 












5 


0.885 to 0.890 


2 


28 


5 


2 


4 


1 


-14 


-3 


-2 


-2 












0.896 to 0.899- 


2 


29 


5 


2 


3 


1 


-14 


-3 


-1 


-1 


-2 


_2 








0.899" to 0.904 


2 


30 


4 


2 


3 


1 


-14 


-4 


-1 


-1 


-2 










0.904 to 1.000 


2 


31 


4 


1 


3 


1 


-14 


-3 


-1 


-1 


-2 


-2 








0.850 to 0.857 


2 


28 


3 


3 


3 


3 


1 


-11 


-2 





-1 










0.857 toO. 809 


2 


28 


3 


3 


4 


2 


1 


-11 


-3 








_1 








0.869 to 0.885 


2 


28 


•1 


3 


4 


2 





-11 


-4 








_1 









0.885 to 0.896 


2 


28 


5 


3 


4 


1 





-10 


-4 


-2 





_1 








0.896 to 0.899- 


2 


29 


5 


3 


3 


L 





-10 


-3 


-2 


-2 


— 1 








0.899" to 0.904 


2 


30 


1 


•1 


2 


1 





-11 


-3 


— 2 


2 










0.904 to 1.000 





31 


i 


3 


2 


1 





-10 


-3 


-2 


-2 


-2 








0.850 to 0.857 


3 


28 


-1 


3 


3 


4 


2 


1 


3 


3 


2 


2 








0.857 to 0.869 


4 


29 


4 


3 


4 


2 


2 


1 


2 


3 


2 


1 








0.869 to 0.885 


4 


29 


5 


3 


4 


2 


1 


1 


1 


3 


2 


1 






11 


0.885 to 0.896 


4 


29 





3 


4 


1 


1 


2 


1 


1 


2 


1 








0.896 to 0.899- 


4 


30 


6 


3 


3 


1 


1 


2 


2 


1 





1 








0.899" to 0.904 


4 


31 


5 


4 


2 


1 


1 


2 


2 


1 





2 








0.904 to 1.000 


4 


32 


5 


3 


2 


1 


1 


2 


2 


1 








12 




0.850 to 0.856 


1 


30 


■> 


-23 


-1 


1 


-I 


-3 


1 


2 


-1 


-1 








0.856 to 0.857 


1 


30 


2 


-23 








-1 


-3 


2 


2 


2 


_2 


-1 






0.857 to 0.869 


1 


30 


2 


-23 


1 


-1 


_2 


-2 


1 


2 





-1 


_1 






O.809 to 0.881 


1 


30 


3 


-23 





-1 


-3 


-2 


1 


2 


1 


-1 


_1 






0.881 to 0.885 


I 


30 


3 


-22 


-1 


-1 


-3 


_2 





2 


2 


-1 


__ 1 






0.885 to 0.896 


1 


30 


1 


-22 


_2 


_2 


-3 








-1 


2 


1 


_1 






0.896 to 0.899- 


1 


31 


4 


-23 


-3 


-2 


-2 





1 


-1 


-1 


1 








0.899" to 0.904 


1 


32 


3 


-23 


-3 


-2 


-2 


-1 


] 





-1 


1 








0.904 to 0.912- 


1 


33 


3 


-24 


-3 


-2 


_2 





1 





-1 


-1 








0.912" to 1.000 


1 


34 


3 


-25 


-3 


-2 


-1 





1 





_2 


-1 








0.850 to 0.856 


1 


30 


3 


2 


-20 





-2 


—2 





1 


1 


-1 








0.856 to 0.857 


1 


30 


3 


2 


-19 


-1 


—2 


-2 


1 


1 





-1 


_1 






0.857 to 0.869 


1 


30 


3 


2 


-18 


-2 


-2 


-2 





1 


1 


-1 


_1 






0.869 to 0.881 


1 


30 


1 


2 


-18 


-3 


-3 


-2 


-1 


2 


1 





_1 




0.SS1 to 0.885 
0.885 to 0.896 


1 


30 


-1 


3 


-19 


-3 


-3 


_2 


-2 


2 


2 





_1 




1 


30 


5 


3 


-19 


-5 


-3 


—1 


-1 





1 


f) 








0.896 to 0.899" 


1 


31 


5 


3 


-20 


-5 


-3 


-1 








-1 











0.899- to 0.904 


1 


32 


1 


1 


-22 


-4 


-3 


-2 











1 








0.904 to 0.912" 


1 


33 


1 


3 


-22 


-4 


-3 


-1 








n 


-1 








0.912- to 1.000 


1 


34 


■1 


2 


-22 


-4 


-2 


-1 








-l 


-1 


-1 
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Table IV C — Continued 





















Ci*(m,») 












n 


m 


g-interval 


X 





























— 






1 


9 


? 2 


?' 


g* 


q> 


<fi 


<?' 


q» 


?» 


gI0 


7" 


12 


0.850 to 0.856 


2 


31 


3 


2 


2 


-16 


-3 


-3 

















0.856 to 0.857 


2 


31 


3 


2 


3 


-17 


-3 


-3 


1 




-1 





-1 






0.857 to 0.869 


2 


31 


3 


2 


4 


-18 


-3 


-3 













-1 






0.860 to 0.881 


2 


31 


4 


2 


4 


-18 


-5 


-3 


-1 




1 








4 


0.881 to 0.885 


2 


31 


4 


3 


3 


-18 


-5 


-3 


-2 




2 








0.885 to 0.896 


2 


31 


5 


3 


3 


-19 


-5 


-2 


-2 


-1 


2 










0.896 to 0.899" 


2 


32 


5 


3 


2 


-19 


-5 


-2 


-1 


-1 












0.899" to 0.904 


2 


33 


4 


4 


1 


-19 


-5 


-3 


-1 


— 1 













0.904 to 0.912" 


2 


34 


4 


3 


1 


-19 


-5 


-2 


-1 


—1 





—1 








0.912 - to 1.000 


2 


35 


4 


2 


1 


-19 


-4 


-2 


-1 


—1 


-1 


_1 


-1 






0.850 to 0.856 


2 


31 


3 


3 


3 


4 


1 


-14 


-2 





-1 


-1 


1 






0.856 to 0.857 


2 


31 


3 


3 


4 


3 


1 


-14 


-1 





-2 


—1 








0.857 to 0.869 


2 


31 


3 


3 


5 


2 


1 


-14 


-2 





-1 


—1 








0.869 to 0.881 


2 


31 


4 


3 


5 


2 





-14 


-3 





-1 


-1 








0.881 to 0.885 


2 


31 


4 


4 


4 


2 





-14 


-4 








-1 






6 


0.885 to 0.896 


2 


31 


5 


4 


4 







-13 


-4 


-2 





-1 








0.896 to 0.899 - 


2 


32 


5 


4 


3 







-13 


-3 


-2 


-2 


—1 








0.899 - to 0.904 


2 


33 


4 


5 


2 







-14 


-3 


-2 


-2 










0.904 to 0.912- 


2 


34 


4 


4 


2 







-13 


-3 


-2 


-2 


-2 








0.912" to 1.000 


2 


35 


4 


3 


2 




1 


-13 


-3 


-2 


-3 


-2 


-1 






0.850 to 0.856 


2 


31 


3 


3 


3 


4 


2 


1 


-9 


-2 


-1 


-1 








0.856 to 0.857 


2 


31 


3 


3 


4 


3 


2 


1 


-8 


-2 


-2 


—1 


—1 






0.857 to 0.869 


3 


31 


4 


3 


5 


2 


1 





-9 


-2 


-1 


—1 


— 1 






0.869 to 0.881 


3 


31 


5 


3 


5 


2 








-10 


-2 


-1 


-1 


— 1 






0.881 to 0.885 


3 


31 


5 


4 


4 


2 








-11 


-2 





—1 


— 1 




7 


0.885 to 0.896 


3 


31 


6 


4 


4 


1 





1 


-11 


-4 





—1 


— 1 






0.896 to 0.899" 


3 


32 


6 


4 


3 


1 





1 


-10 


-4 


-2 


—1 


— 1 






0.899" to 0.904 


3 


33 


5 


5 


2 


1 








-10 


-4 


-2 





— 1 






0.904 to 0.912- 


3 


34 


5 


4 


2 


1 





1 


-10 


-4 


-2 


-2 


— 1 






0.912- to 1.000 


3 


35 


5 


3 


2 


1 


1 


1 


-10 


-4 


-3 


-2 


-2 






0.850 to 0.856 


4 


32 


4 


3 


3 


4 


2 


1 


3 


3 


2 


2 


3 






0.856 to 0.857 


4 


32 


4 


3 


4 


3 


2 


1 


4 


3 


1 


2 


2 






0.857 to 0.869 


4 


32 


4 


3 


5 


2 


2 


1 


3 


3 


2 


2 


2 






0.869 to 0.881 


4 


32 


5 


3 


5 


2 




1 


2 


3 


2 


2 


2 






0.881 to 0.885 


4 


32 


5 


4 


4 


2 




1 


1 


3 


3 


2 


2 




12 


0.885 to 0.896 


4 


32 


6 


4 


4 


1 




2 


1 




3 


2 


2 






0.896 to 0.899- 


4 


33 


6 


4 


3 


1 




2 


2 




1 


2 


2 






0.899" to 0.904 


4 


34 


5 


5 


2 


2 




1 


2 




1 


2 


1 






0.904 to 0.912- 


4 


35 


5 


4 


2 


2 




2 


2 




1 





1 






0.912- to 1.000 


4 


36 


5 


3 


2 


2 


2 


2 


2 











The exponential symbols +, - indicate only the relative magnitude of two 
different roots that are equal to three decimal places (i.e., a < a + ). 



Table V A — Expected Number of Tests Required and Size x of 

the Next Sample To Be Taken for Any Gi(m,n) and Any 

Hi(n) Situation When Procedure #1 is Used and q = 0.90 





AM 




Gi(w,b) 


It 


X 
























m = 2 


»i = 3 


m = 4 


IB = 5 


m = 6 


m = 7 


m = 8 


m = 9 








* = 1 


x = 1 


x = 2 


1 = 2 


x = 2 


x = 3 


x = 4 


x = 4 


1 


1.000 


1 


















2 


1.290 


2 


1.526 
















3 


1.661 


3 


2.153 


2.439 














4 


2.051 


4 


2.485 


2.971 


3.056 












5 


2.490 


5 


2.866 


3.325 


3.547 


3.637 










6 


2.943 


6 


3.282 


3.727 


3.925 


4.115 


4.148 








7 


3.414 


7 


3.728 


4.157 


4.343 


4.512 


4.621 


4.628 






8 


3.904 


8 


4.191 


4.613 


4.784 


4.943 


5.033 


5.100 


5.100 




!) 


4.395 


9 


4.672 


5.085 


5.250 


5.396 


5.476 


5.528 


5.575 


5.543 


10 


4.872 


6 


5.163 


5.570 


5.728 


5.868 


5.938 


5.980 


6.014 


6.021 


11 


5.327 


6 


5.646 


6.056 


6.210 


6.346 


6.411 


6.445 


6.470 


6.466 


12 


5.790 


6 


6.111 


6.528 


6.687 


6.821 


6.883 


6.914 


6.933 


6.922 


13 


6.261 


7 


6.570 


6.993 


7.157 


7.294 


7.356 


7.385 


7.401 


7.385 


14 


6.732 


7 


7.037 


7.456 


7.623 


7.765 


7.829 


7.857 


7.872 


7.853 


15 


7.213 


7 


7.509 


7.925 


8.089 


8.232 


8.300 


8.330 


8.345 


8.324 


16 


7.695 


7 


7.985 


8.400 


8.561 


8.702 


8.770 


8.802 


8.819 


8.798 


17 


8.161 


6 


8.467 


8.878 


9.038 


9.176 


9.242 


9.275 


9.293 


9.273 


18 


8.629 


6 


8.940 


9.354 


9.513 


9.651 


9.715 


9.746 


9.764 


9.746 


19 


9.100 


7 


9.407 


9.825 


9.986 


10.124 


10.188 


10.218 


10.235 


10.217 


20 


9.572 


7 


9.877 


10.294 


10.458 


10.596 


10.660 


10.000 


10.706 


10.687 


21 


10.044 


7 


10.348 


10.764 


10.927 


11.068 


11.133 


11.163 


11.179 


11.159 


22 


10.520 


7 


10.820 


11.236 


11.398 


11.538 


11.605 


11.635 


11.651 


11.632 


23 


10.996 


6 


11.295 


11.709 


11.871 


12.010 


12.076 


12.108 


12.124 


12.105 


24 


11.466 


6 


11.770 


12.184 


12.345 


12.484 


12.549 


12.580 


12.597 


12.578 


25 


11.937 


7 


12.243 


12.658 


12.819 


12.957 


13.022 


13.052 


13.069 


13.051 


26 


12.408 


7 


12.714 


13.130 


13.292 


13.430 


13.495 


13.525 


13.541 


13.522 


27 


12.881 


7 


13.185 


13.601 


13.763 


13.902 


13.967 


13.997 


14.014 


13.994 


2S 


13. 353 


7 


13.657 


14.073 


14.235 


14.374 


14.440 


14.470 


14.486 


14.467 


2!) 


13.827 


7 


14.129 


14.545 


14.707 


14.846 


14.912 


14.943 


14.959 


14.939 


30 


14.301 


7 


14.603 


15.018 


15.179 


15.319 


15.384 


15.415 


15.431 


15.412 


31 


14.773 


7 


15.077 


15.491 


15.653 


15.792 


15.857 


15.888 


15.904 


15.885 


32 


15.244 


7 


15.549 


15.965 


16.126 


16.264 


16.329 


16.360 


16.377 


16.358 


33 


15.717 


7 


16.021 


16.437 


16.598 


16.737 


16.802 


16.833 


16.849 


16.830 


34 


16.189 


7 


16.493 


16.909 


17.071 


17.210 


17.275 


17.305 


17.321 


17.302 


35 


16.661 


7 


16.965 


17.381 


17.543 


17.682 


17.747 


17.778 


17.794 


17.775 


36 


17.135 


7 


17.438 


17.853 


18.015 


18.154 


18.219 


18.250 


18.266 


18.247 


37 


17.608 


7 


17.910 


18.326 


18.487 


18.626 


18.692 


18.723 


18.739 


18.720 


38 


18.080 


7 


18.384 


18.799 


18.960 


19.099 


19.164 


19.195 


19.212 


19.192 


30 


18.552 


7 


18.856 


19.272 


19.433 


19.572 


19.637 


19.668 


19.684 


19.665 


40 


19.024 


7 


19.328 


19.744 


19.906 


20.044 


20.110 


20.140 


20.157 


20.137 


41 


19.497 


7 


19.801 


20.216 


20.378 


20.517 


20.582 


20.613 


20.629 


20.610 


42 


19.969 


7 


20.273 


20.68S 


20.850 


20.989 


21.055 


21.085 


21 . 102 


21.082 


43 


20.442 


7 


20.745 


21.161 


21.323 


21.462 


21.527 


21.558 


21.574 


21.555 


44 


20.915 


7 


21.218 


21.634 


21.795 


21.934 


21.999 


22.030 


22.047 


22.027 


45 


21.387 


7 


21.691 


22.106 


22.268 


22.407 


22.472 


22.503 


22.519 


22.500 


46 


21.860 


7 


22.164 


22.579 


22.740 


22.879 


22.945 


22.975 


22.992 


22.972 


47 


22.332 


7 


22.636 


23.051 


23.213 


23.352 


23.417 


23.448 


23.464 


23.445 


48 


22.805 


7 


23 . 108 


23.524 


23.685 


23.825 


23.890 


23.920 


23.937 


23.917 


49 


23.277 


7 


23.581 


23.996 


24.158 24.297 


24.362 


24.393 


24.409 


24.390 


50 


23.750 


7 


24.053 


24.469 


24.630J 24.769 


24.835 


24.865 


24.882 


24.862 
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Table V A — Continued 





ffl(ll) 


X 


Gi(m,n) 


n 


m = 2 


m = 3 


HI = 4 


m = 5 


m = 6 


m = 7 


in = 8 


M -= 9 








x = 1 


x = I 


x = 2 


x = 2 


* = 2 


x = 3 


* = 4 


1 = 4 


51 


24.222 


7 


24.526 


24.941 


25.103 


25.242 


25.307 


25.338 


25.354 


25.335 


52 


24 . 695 


7 


24.998 


25.414 


25.575 


25.714 


25.780 


25.810 


25.827 


25.807 


53 


25.167 


7 


25.471 


25.886 


26.048 


26.187 


26.252 


26.283 


26.299 


26.280 


54 


25.640 


7 


25.943 


26.359 


26.521 


26.660 


26.725 


26.755 


26.772 


26.752 


55 


26.112 


7 


26.416 


26.831 


26.993 


27.132 


27.197 


27.228 


27.244 


27.225 


56 


26.585 


7 


26.888 


27.304 


27.465 


27.605 


27.670 


27.700 


27.717 


27.697 


57 


27.057 


7 


27.361 


27.776 


27.938 


28.077 


28.142 


28.173 


28.189 


28.170 


58 


27.530 


7 


27.833 


28.249 


28.410 


28.550 


28.615 


28.645 


28.662 


28.643 


59 


28.002 


7 


28.306 


28.721 


28.883 


29.022 


29.087 


29.118 


29.134 


29.115 


60 


28.475 


7 


28.779 


29.194 


29.356 


29.495 


29.560 


29.590 


29.607 


29.588 


61 


28.947 


7 


29.251 


29.666 


29.828 


29.967 


30.032 


30.063 


30.079 


30.060 


62 


29.420 


7 


29.723 


30.139 


30.301 


30.440 


30.505 


30.535 


30.552 


30.533 


63 


29.892 


7 


30.196 


30.611 


30.773 


30.912 


30.977 


31.008 


31.024 


31.005 


64 


30.365 


7 


30.668 


31.084 


31.246 


31.385 


31.450 


31.480 


31.497 


31.478 


65 


30.837 


7 


31.141 


31.556 


31.718 


31.857 


31.922 


31.953 


31.969 


31.950 


66 


31.310 


7 


31.614 


32.029 


32.191 


32.330 


32.395 


32.425 


32.442 


32.423 


67 


31.782 


7 


32.086 


32.502 


32.663 


32.802 


32.867 


32.898 


32.914 


32.895 


68 


32.255 


7 


32.559 


32.974 


33.136 


33.275 


33.340 


33.370 


33.387 


33.368 


69 


32.727 


7 


33.031 


33.446 


33.608 


33.747 


33.812 


33.843 


33.859 


33.840 


70 


33.200 


7 


33.504 


33.919 


34.081 


34.220 


34.285 


34.316 


34.332 


34.313 


71 


33.672 


7 


33.976 


34.392 


34.553 


34.692 


34.757 


34.788 


34.804 


34.785 


72 


34.145 


7 


34.449 


34.864 


35.026 


35.165 


35.230 


35.261 


35.277 


35.258 


73 


34.617 


7 


34.921 


35.337 


35.498 


35.637 


35.702 


35.733 


35.749 


35.730 


74 


35.0iK) 


7 


35.394 


35.809 


35.971 


36.110 


36.175 


36.206 


36.222 


36.203 


75 


35.562 


7 


35.866 


36.282 


36.443 


36.582 


36.647 


36.678 


36.694 


36.675 


76 


36.035 


7 


36.339 


36.754 


36.916 


37.055 


37.120 


37.151 


37.167 


37.148 


77 


36.507 


7 


36.811 


37.227 


37.388 


37.527 


37.592 


37.623 


37.639 


37.620 


78 


36.980 


7 


37.284 


37.699 


37.861 


38.000 


38.065 


38.096 


38.112 


38.093 


79 


37.453 


7 


37.756 


38.172 


38.333 


38.472 


38.537 


38.568 


38.584 


38.565 


80 


37.925 


7 


38.229 


38.644 


38.806 


38.945 


39.010 


39.041 


39.057 


39.038 


81 


38.398 


7 


38.701 


39.117 


39.278 


39.417 


39.482 


39.513 


39.530 


39.510 


82 


38.870 


7 


39.174 


39.589 


39.751 


39.890 


39.955 


39.986 


40.002 


39.083 


83 


39.343 


7 


39.646 


40.062 


40.223 


40.362 


40.427 


40.458 


40.475 


40.455 


84 


39.815 


7 


40.119 


40.534 


40.696 


40.835 


40.900 


40.931 


40.947 


40.928 


85 


40.288 


7 


40.591 


41.007 


41.168 


41.307 


41.373 


41.403 


41.420 


41.400 


86 


40.760 


7 


41.064 


41.479 


41.641 


41.780 


41.845 


41.876 


41.892 


41.873 


87 


41.233 


7 


41.536 


41.952 


42.113 


42.252 


42.318 


42.348 


42.365 


42.345 


ss 


41.705 


7 


42.009 


42.424 


42.586 


42.725 


42.790 


42.821 


42.837 


42.818 


89 


42.178 


7 


42.481 


42.897 


43.058 


43.197 


43.263 


43.293 


43.310 


43.290 


90 


42.650 


7 


42.954 


43.369 


43.531 


43.670 


43.735 


43.766 


43.782 


43.763 


•H 


43.123 


7 


43.426 


43.842 


44.003 


44.142 


44.208 


44.238 


44.255 


44.235 


92 


43.595 


7 


43.899 


44.314 


44.476 


44.615 


44.680 


44.711 


44.727 


44.708 


93 


44.068 




44.371 


44.787 


44.948 


45.087 


45.153 


45.183 


45.200 


45.180 


94 


44.540 


7 


44.844 


45.259 


45.421 


45.560 


45.625 


45.656 


45.672 


45.653 


95 


45.013 




45.316 


45.732 


45.893 


46.032 


46.098 


46.128 


46.145 


46.125 


96 


45.485 


< 


45.789 


46.204 


46.366 


46.505 


46.570 


46.601 


46.617 


46.598 


( .)7 


45.958 




46.261 


46.677 


46.838 


46.978 


47.043 


47.073 


47.090 


47.070 


98 


46.43C 


7 


46.734 


47.14S 


47.311 


47.450 


47.515 


47.546 


47.562 


47.543 


'.ID 


46.903 


7 


47.206 


47.622 


47.784 


47.923 


47.988 


48.018 


48.035 


48.016 


100 


47.375 


7 


47.675 


48.09-1 


48.256 


48.395 


48.46C 


48.491 


48.507 


48.488 
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Table VI — Expected Number of Tests Required for Procedure R 2 
The integer below q" opposite H 2 (n) is the coefficient of q" in the polynomial 
formula for H 2 (n). 





g-interval 


-V 


l 

2 
3 


Q 
-1 


-1 


„3 


g* 


</ 6 


<7 6 


s 7 


9 8 


1 

-1 

2 

-1 

-1 

-1 

4 

-1 


glO 


9" 


,U 


ff*(2) 


0.000 to 0.618 
0.618 to l.O(M) 


1 

2 




ff*(3) 


0.000 to 0.618 
0.618 to 0.755 
0.755 to 1.000 


1 
2 
3 


3 
5 
5 


-3 
-2 


-1 
-1 


1 
-1 














ff 3 (4) 


0.000 to 0.618 
0.618 to 0.755 
0.755 to 0.819 
0.819 to 1.000 


1 
2 
3 

4 


4 

7 
7 
8 


-5 
-3 

-4 



-2 
-2 


1 
-2 
-1 


-1 
2 







— 






ff*(5) 


0.000 to 0.618 
0.618 to 0.755 
0.755 to 0.819 
0.819 to 0.857 
0.857 to 1.000 


1 
2 
•'5 
4 
5 


5 

9 

9 

11 

11 


-7 
-4 
-7 
-7 


1 

-3 
-2 
-2 



-2 

-1 



-1 
3 

-1 



1 

-1 

2 

-1 




ff*(6) 


0.000 to 0.618 
0.618 to 0.755 
0.755 toO.Sli) 
0.819 to 0.857 
0.857 to 0.881 
0.881 to 0.885 
0.885 to 1.000 


1 
2 
3 
4 
5 
6 
6 


6 
11 
11 
14 
14 
14 
14 


-9 

-5 

-10 

-10 

-10 

-9 


2 

-4 
-2 
-2 
-1 
-2 


-1 

-2 
-1 

1 



-1 




4 

-1 

-1 

-1 




1 
-1 

3 
-3 
-1 
-1 


-2 

4 
-4 
-1 
-2 
-2 


-1 
-1 
-1 

3 

-1 
-2 
-1 
5 
-2 
-2 


-2 
-2 

7 
-3 
-3 



1 

1 
-2 
-2 

8 

-4 

-4 



2 

1 



1 

2 

-1 

-1 

3 

3 

1 

3 

-1 

-1 

5 

5 

-2 





4 

-1 

-2 

7 
7 
-3 
-1 

1 


-1 

-1 

-1 
-1 
-1 

3 
-1 

-1 
-1 
-1 
-1 

-1 
-1 

5 
-4 
-1 






H,(7) 


0.000 to 0.618 
0.618 to 0.755 
0.755 to 0.819 
0.819 to 0.857 
0.857 to 0.881 
0.881 to 0.885 
0.885 to 0.899- 
0.899" to 1.000 


1 

2 
3 
4 
5 
6 
6 
7 


7 
13 
13 
17 
17 
17 
17 
17 


-11 
-6 
-13 
-13 
-13 
-11 
-11 


3 
-5 
-2 
-2 


-3 
-3 


-2 

-2 

-1 

2 

-1 

-2 

-1 


1 

5 
-1 
-2 

-1 
2 

1 




ff.(8) 


0.000 to 0.618 
0.618 to 0.755 
0.755 to 0.819 
0.819 to 9.857 
0.857 to 0.881 
0.881 to 0.885 
0.885 to 0.899- 
0.899- to 0.912 
0.912 to 0.932 
0.932 to 1.000 


1 

2 
3 
4 
5 
6 
6 
7 
8 
8 


8 
15 
15 
20 
20 
20 
20 
20 
20 
20 


-13 

-7 
-16 
-16 
-16 
-13 
-13 
-13 
-14 


4 
-6 
-2 
-2 

1 
-4 
-4 
-4 
-4 


-3 
-2 

-1 

3 

-2 

-3 

-1 
-1 
-1 


2 
6 


-3 

-1 

4 
1 
1 

1 


-1 
-3 
4 
-5 
-1 
-3 
-2 
-2 
-2 




ff 2 (9) 


0.000 to 0.618 
0.618 to 0.755 
0.755 to 0.819 
0.819 to 0.857 
0.857 to 0.881 
0.SS1 to().S9() 
0.S90 to 0.899- 
0.899- to 0.912 
0.912 to 0.922 
0.922 to 0.938 
0.938 to 1.000 


1 

2 
3 

4 
5 

li 

r> 

7 
8 
9 
9 


9 
17 
17 
23 
23 
23 
23 
23 
23 
24 
25 


-15 

-8 
-19 
-19 
-19 
-15 
-15 
-15 
-16 
-18 


5 
-7 
—2 
-2 

2 

-5 
-5 
-5 
-5 

-4 


-4 
-2 
-1 

4 
-3 

-4 

-1 
-1 
-1 

-1 


3 

7 
1 

-4 

-1 
6 
1 
1 
1 
1 


-2 
-4 

3 
-5 
-1 

-4 
-2 
-2 
-2 

-2 
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Table 


VI 


— Continued 


















Q- 


nterval 


X 

1 


l 
10 


Q 


Q*- 


:/ 3 


Q A 


9 s 


q* 


V 7 


'i' 


<7» 


,10 


q" 


</'•' 




0.000 


to 0.618 






0.618 


to 0.755 


2 


19 


-17 


6 


-5 


4 


-3 


2 


-1 





1 


-1 








0.755 


to 0.819 


3 


19 


-9 


-8 


-2 


8 


-5 


-2 


5 


-2 


-2 


2 








0.819 


to 0.857 


4 


26 


— 22 


-2 


-1 


2 


2 


-2 


-1 


-1 


3 


-1 








0.857 


to 0.881 


5 


26 


-22 


-2 


5 


-5 


-5 


8 


-2 








-1 








0.881 


to 0.890 


(J 


26 


-22 


3 


-1 


-1 


-1 


-4 


8 


-2 


-1 








ffi(10) 


0.890 


to 0.899- 


6 


26 


-17 


-6 


-5 


8 


-5 


-4 


8 


-2 


-1 










0.899- 


to 0.912 


7 


26 


-17 


-6 


-1 




-2 





-4 


7 


-1 


-1 








0.912 


to 0.922 


8 


26 


-17 


-6 


-1 




-2 


3 


-2 


-(i 


7 


-1 








0.922 


to 0.930 


9 


28 


-20 


-5 


-1 




-2 


1 


1 


-2 


-3 


4 








0.930 


to 0.938 


10 


28 


-20 


-5 


-1 




-1 


1 


-1 


-1 












0.938 


to 0.900 


10 


29 


-22 


-4 


-1 




-1 











-1 










0.960 


to 1.000 


10 


29 


-22 


-4 







-3 





1 





-1 










0.000 


to 0.018 


1 


11 




























0.618 


to 0.755 


2 


21 


-19 


7 


-6 


5 


-4 


3 


-2 


1 





-1 


1 






0.755 


to 0.819 


3 


21 


-10 


-9 


-2 


9 


-6 


-2 


6 


-3 


-2 


3 


-1 






0.819 


to 0.857 


4 


29 


-25 


-2 


-1 


3 


1 


-2 


-1 


-1 


4 


-1 


-1 






0.857 


to 0.881 


5 


29 


-25 


-2 


6 


-6 


-5 


8 


-2 


1 


-1 


-3 


3 






0.881 


to 0.890 


6 


29 


-25 


4 


-5 


-1 


-1 


-4 


8 


-2 








-1 




//,(11) 


0.890 


to 0.899- 


6 


29 


-19 


-7 


-6 


10 


-0 


-4 


8 


-2 








-1 




0.899" 


to 0.912 


7 


29 


-19 


-7 


-1 




-2 





-4 


8 


-2 


-1 








0.912 


to 0.922 


8 


29 


-19 


-7 


-1 




-2 


4 


-3 


-8 


10 


-1 


-1 






0.922 


to 0.930 


9 


32 


-24 


-5 


-1 




-2 


1 


2 


-3 


-5 


7 


-1 






0.930 


to 0.930 


10 


32 


-24 


-5 


-1 










-2 


1 


-1 


3 


■1 






0.936 


to 0.938 


11 


32 


-24 


-5 







-2 


1 


-1 


-1 












0.938 


to 0.900 


11 


33 


-26 


-4 







-2 











-1 










0.960 


to 1.000 


11 


33 


-26 


-4 


1 





-4 


2 


1 


-1 


-1 










0.000 


to 0.018 


1 


12 




























0.018 


to 0.755 


2 


23 


-21 


8 


— 7 


6 


-5 


4 


-3 


2 


-1 





1 


-1 




0.755 


to 0.819 


3 


23 


-11 


-10 


-2 


10 


—7 


-2 


7 


-4 


-2 


4 


-1 


-1 




0.819 


to 0.857 


4 


32 


-28 


-2 


-1 


4 





-2 


-1 





4 


-2 


-1 






0.S57 


to 0.881 


5 


32 


-28 


-2 


7 


-7 


-5 


8 


-2 


2 


-2 


-4 


5 


-1 




0.SS1 


to 0.890 


6 


32 


-28 


5 


-6 


-1 


-1 


-4 


8 


-1 





-1 


-1 






0.890 


to 0.899- 


6 


32 


-21 


-8 


-7 


12 


-7 


-4 


9 


-2 


-1 





-1 




// 2 (12) 


0.899" 


to 0.912 


7 


32 


-21 


-8 


-1 


1 


-2 





-4 


8 


-2 








-1 


0.912 


to 0.922 


8 


32 


-21 


-8 


-1 


1 


-2 


5 


-4 


-9 


12 


-2 


-1 






0.922 


to 0.930 


9 


36 


-28 


-5 


-1 


1 


_2 


1 


3 


-4 


-7 


10 


-1 


-1 




0.930 


to 0.930 


10 


30 


-2S 


-5 


-1 


1 


1 


-1 


-3 


3 


-2 


-5 


7 


-1 




0.930 


to 0.938 


11 


30 


-28 


— 5 


1 





-3 


3 


-2 


-2 


2 


-1 


-3 


4 




0.93S 


to 0.941 


11 


37 


-30 


-4 


1 





-3 


2 


-1 


-1 


1 


-1 


-3 


4 




0.941 


to 0.900 


12 


37 


-30 


-4 


1 





-3 


2 


-1 


I) 





-1 








0.900 


to 0.972 


12 


37 


-30 


-4 


2 


-1 


-5 


4 


o 


-1 





-1 








0.972 


to 1.000 


12 


37 


-29 


-4 








-4 


2 


1 


-1 


-1 









The exponential symbols +, — indicate only the relative magnitude of two 
different roots that are equal to three decimal places (i.e., a~ < a + ). 



Table VII — The Dividing Points Between x and x + 1 for the 
Information Procedure R-> 

G-Situation 











m 


























4 


6 


7 


8 


9 


11 


12 


13 


14 


16 


1 


0.7549 


0.6518 


0.6369 


0.6289 


0.6245- 


0.6204 


0.6195" 


0.6189 


0.6186 


0.6182 






0.8899 


0.8376 


0.8087 


0.7913 


0.7728 


0.7077 


0.7642 


0.7617 


0.7580 










0.9378 


0.9016 


0.8631 


0.8524 


0.8446 


0.8388 


0.8313 














0.9340 


0.9160 


0.0031 


0.8935 


0.8806 


5 
6 
7 














0.9723 


0.9528 


0.9384 


0.9193 


















0.9796 


0.0529 




















0.9844 



(j-Situation 



hi 




10 


20 


30 


40 


50 


60 


70 


80 


90 


100 




0.6219 


0.6181 


0.6180 


0.6180 


0.6180 


0.6180 


0.6180 


0.6180 


0.6180 


0.6180 


1 


0.7802 


0.7560 


0.7549 


0.75 19 


0.7549 


0.7549 


0.7549 


0.7549 


0.7549 


0.7549 


2 


0.8786 


0.8241 


0.8198 


0.8193 


0.8192 


0.8192 


0.8192 


0.8192 


0.8192 


0.8192 


3 


0.9601 


0.8677 


0.8587 


0.8571 


0.8568 


0.8567 


0.8567 


0.8567 


0.8567 


0.8567 


1 




0.8999 


0.8854 


0.8823 


0.8816 


0.8814 


0.8813 


0.8813 


0.8813 


0.8813 


5 




0.9262 


0.9056 


0.9008 


0.8993 


0.8989 


0.8987 


0.8987 


0.8987 


0.8987 







0.9490 


0.9218 


0.9150 


0.9129 


0.9121 


0.9118 


0.9117 


0.9116 


0.9116 


7 




0.9699 


0.9355" 


0.9267 


0.9236 


0.9225- 


0.9220 


0.9218 


0.9217 


0.9216 


8 




0.9900 


0.9474 


0.9364 


0.9325+ 


0.9309 


0.9302 


0.9299 


0.9297 


0.9296 


9 






0.9582 


0.9449 


0.9400 


0.9380 


0.9370 


0.9365* 


0.9303 


0.0362 


10 






0.9682 


0.9524 


0.9466 


0.9440 


0.9428 


0.9422 


0.9418 


0.9417 


11 






0.9776 


0.9593 


0.9524 


0.9493 


0.9478 


0.9470 


0.9466 


0.9403 


Pi 






0.9867 


0.9655+ 


0.9576 


0.9540 


0.9522 


0.9512 


0.9507 


0.9504 


13 




0.9991 


0.9956 


0.9714 
0.9770 


0.9624 
0.9668 


0.9582 
0.9620 


0.9561 
0.9590 


0.9549 
0.9583 


0.9543 
0.9575" 


0.9539 
0.9570 


14 


46 


15 


45 


0.99S2 




0.9823 


0.9709 


0.9656 


0.9028 


0.9613 


0.9604 


0.9598 


16 


44 


0.9973 




0.9874 


0.9747 


0.9689 


0.9658 


0.9010 


0.9630 


0.9023 


17 


43 


0.9961 




0.9925- 


0.9784 


0.9719 


0.9685- 


0.9665+ 


0.9654 


0.9610 


18 


42 


0.9955' 




0.9975" 


0.9820 


0.9748 


0.9710 


0.9089 


0.9676 


0.0607 


10 


41 


0.9946 


0.9989 




0.9854 


0.9776 


0.9734 


0.9711 


0.9090 


0.9687 


20 


40 


0.9937 


0.9978 




0.9887 


0.9802 


0.0757 


0.9731 


0.9715" 


0.9705- 


21 


30 


0.9928 


0.9967 




0.9920 


0.9827 


0.9779 


0.9750* 


0.9733 


0.9721 


22 


38 


0.9919 


0.9956 




0.9952 


0.9852 


0.9799 


0.9769 


0.9750" 


0.9737 


23 


37 


0.9909 


0.9944 




0.9984 


0.9876 


0.9819 


0.9780 


0.9705' 


0.9752 


24 


36 


0.9899 


0.9933 


0.9986 




0.9899 


0.9838 


0.9803 


0.9780 


0.9766 


25 


35 


0.9889 


0.9921 


0.9972 




0.9922 


0.9856 


0.9819 


0.9795- 


0.9779 


20 


34 


0.9879 


0.9910 


0.9957 




0.9944 


0.9874 


0.9834 


0.9S09 


0.9702 


27 


33 


0.9869 


0,9898 


0.9943 




0.9967 


0.9892 


0.9849 


0.9822 


0.9804 


28 


32 


0.9858 


0.9886 


0.9928 




0.9989 


0.9909 


0.9863 


0.9835- 


0.9816 


29 


31 


0.9847 


0.9873 


0.9914 


0.9981 




0.9926 


0.9877 


0.9847 


0.9827 


30 


30 


0.9836 


0.9860 


0.9899 


0.9962 




0.9943 


0.9891 


0.9859 


0.9838 


:u 


29 


0.9824 


0.9847 


0.9883 


0.9943 




0.9959 


0.9904 


0.9871 


0.9848 


32 


28 


0.9812 


0.9834 


0.9868 


0.9924 




0.9976 


0.9918 


0.9882 


0.9858 


33 


27 


0.9800 


0.9820 


11.9852 


0.9904 




0.9992 


0.9931 


0.9893 


0.9868 


34 


26 


0.9786 


0.9805+ 


0.9835+ 


0.9884 


0.9974 




0.9943 


0.9904 


0.9878 


35 


25 


0.9773 


0.9790 


0.9818 


0.9864 


0.9947 




0.9956 


0.9914 


0.9887 


36 


24 


0.9758 


0.9775- 


O.980O 


0.9843 


0.9920 




0.9969 


0.9925" 


0.9896 


37 


23 


0.9743 


0.9758 


0.9782 


0.9822 


0.9893 




0.9981 


0.9935+ 


0.9905*- 


38 


22 


0.9727 


0.9741 


0.9763 


0.9800 


0.9866 




0.9994 


0.9945* 


0.9914 


39 


21 


0.9709 


0.9722 


0.9742 


0.9776 


0.9837 


0.9961 




0.9955+ 


0.9922 


40 


20 


0.9691 


0.9702 


0.9721 


0.9752 


0.9808 


0.9921 




0.9965+ 


0.9931 


41 


19 


0.9671 


0.9681 


0.9698 


0.9727 


0.9778 


0.9881 




0.9975' 


0.9939 


42 


18 


0.964!) 


0.9659 


0.9674 


0.9700 


0.9746 


0.9840 




0.9985+ 


0.9948 


43 


17 


0.9626 


0.9634 


0.9648 


0.9671 


0.9713 


0.979S 




0.9995+ 


0.9950 


44 


16 


0.9600 


0.9608 


0.9619 


0.9640 


0.9678 


0.9754 


0.9935" 




0.9964 


15 


15 


0.9572 


0.9578 


0.9588 


0.9606 


0.9640 


0.9708 


0.9869 




0.9972 


40 


11 


0.9541 


0.9546 


0.9554 


0.9570 


0.9599 


0.9659 


0.9802 




0.9980 


47 


13 


0.9505" 


0.9509 


0.9516 


0.9529 


U.9555 


0.9607 


0.9733 




0.998S 


48 


12 


0.9464 


0.9467 


0.9473 


0.9484 


0.9506 


0.9551 
0.9489 


0.9660 
0.9582 


0.9872 


0.9990 


49 


11 


0.9417 


0.9420 


0.9424 


0.9433 


0.9451 






10 


0.9362 


0.9364 


0.9367 


0.9374 


0.9388 


0.9419 


0.9498 


0.9742 






9 


0.9297 


0.9298 


0.9300 


0.9305" 


0.9315+ 


0.9340 


0.9404 


0.960S 






S 


0.9216 


0.9217 


0.9218 


0.9222 


0.9229 


0.9248 


0.9298 


0.9463 






7 


0.9116 


0.9116 


0.9117 


0.9119 


0.9124 


0.9137 


0.9174 


0.9303 






B 


0.8987 


0.8987 


0.8987 


0.8988 


0.8991 


0.8999 


0.9024 


0.9120 


0.9645" 




5 


0.8813 


0.8813 


0.8813 


0.8813 


0.8814 


0.8818 


0.8834 


0.8898 


0.9270 




4 


0.8567 


0.8567 


0.8567 


0.8567 


0.8567 


0.8569 


0.8576 


0.8612 


0.8862 




3 


0.8192 


0.8192 


0.8192 


0.8192 


0.8192 


0.8192 


0.8194 


0.8209 


0.8346 




■1 


0.7549 


0.7549 


0.7549 


0.7549 


0.7549 


0.7549 


0.7549 


0.7552 


0.7599 




1 


0.6180 


0.6180 


0.6180 


0.6180 


0.6180 


0.6180 


0.6180 


0.6180 


0.0184 


0.6823 




95 


85 


75 


65 


55 


45 


35 


25 


15 


5 


X 


m 
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Table VII — Continued 
//-Situation 



X 


Q 


X 


Q 


1 


0.6180 


51 


0.9866 


2 


0.7549 


52 


0.9869 


3 


0.8192 


53 


0.9871 


4 


0.8567 


54 


0.9874 


5 


0.8813 


55 


0.9876 


6 


0.8987 


56 


0.9878 


7 


0.9116 


57 


0.9880 


8 


0.9216 


58 


0.9882 


9 


0.9200 


59 


0.9884 


10 


0.9361 


60 


0.9886 


11 


0.0415" 


61 


0.9888 


12 


0.9400 


62 


0.9890 


13 


0.9499 


63 


0.9891 


14 


0.9533 


64 


0.9893 


IS 


0.9563 


65 


0.9895" 


Hi 


0.9588 


66 


0.9896 


17 


0.9612 


67 


0.9898 


IS 


0.9632 


68 


0.9899 


19 


0.9651 


69 


0.9901 


20 


0.9667 


70 


0.9902 


21 


II. WW! 


71 


0.9904 


22 


0.9697 


72 


0.9905" 


23 


0.9709 


73 


0.9906 


24 


0.9721 


74 


0.9907 


25 


0.9732 


75 


0.9909 


26 


0.9742 


76 


0.9910 


27 


0.9751 


77 


0.9911 


28 


0.9701) 


78 


0.9912 


29 


0.9768 


79 


0.9913 


30 


0.9775+ 


80 


0.9914 


31 


0.9782 


81 


0.0915+ 


32 


0.9789 


82 


0. Oil Hi 


33 


0.9795+ 


83 


0.9917 


34 


0.9801 


84 


0.9018 


35 


0.9807 


85 


0.9919 


36 


0.9812 


86 


0.9920 


37 


0.9817 


87 


0.9921 


38 


0.9822 


88 


0.9922 


39 


0.9826 


89 


0.9923 


40 


0.9830 


90 


0.9924 


41 


0.9834 


91 


0.9925" 


42 


0.9838 


92 


0.9925+ 


43 


0.9842 


93 


0.9926 


44 


0.9845+ 


94 


0.9927 


45 


0.9849 


95 


0.9928 


46 


0.9852 


96 


0.9928 


47 


0.9855+ 


97 


0.9929 


48 


0.985S 


98 


0.9930 


49 


0.9801 


99 


0.9931 


50 


0.9864 


100 


0.9931 



The exponents + and ~ are to be used for rounding in the usual manner. 
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Table VIII — Polynomial Expressions for F 2 *(m) 
for Procedure Ri 













Coefficient of 






























j-interval 


Next Test 
Group Size 




i 
1 


9 
1 


9 1 


9 1 


9 4 


9' 


9° 


9' 


9" 


9' 


9'° 


9" 




2 


O.(MK) g g g 1.000 


1 


3 


1 


2 


2 




















O.(MX) < tf < 1.000 


1 


4 


2 


2 


2 


2 


















0.755 ^ (/ < 1.000 


2 


5 


2 


2 


2 


3 


3 
















0.682 < q < 1.000 


2 


•{ 

7 


2 


2 


3 


3 


3 


3 














0.652 < q < 0.890 


2 


2 


3 


3 


2 


3 


3 














0.890 < q ^ 1.000 


3 


2 


3 


3 


3 


3 


3 


3 












0.838 <^ 1.000 


3 


9 


2 


3 


3 


3 


3 


3 


4 


4 










0.809 < q < 0.938 


3 


3 


3 


3 


3 


3 


3 


3 


3 










0.938 < q < 1.000 


4 


3 


3 


3 


3 


3 


3 


3 


4 


4 








0.902 <: q < 1.000 


4 


11 


3 


3 


3 


3 


3 


4 


4 


3 


4 


4 






0.879 < q < 0.960 


4 


3 


3 


3 


4 


4 


3 


3 


3 


4 


4 






0.960 < q < 1-000 


5 


3 


3 


3 


4 


4 


3 


4 


4 


3 


4 


4 




0.934 < q <■ 1.000 


5 


»{ 


3 


3 


3 


4 


4 


3 


4 


4 


4 


4 


4 


4 


0.916 < q ^ 0.972 


5 


3 


4 


4 


3 


4 


4 


3 


4 


4 


3 


4 


4 


0.972 < q < 1.000 


6 



